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I.  INTRODUCTION 


The  1972  Pacific  Marginal  Ice  Zone  Study  (MIZPAC-72)  by  the  Applied 
Physics  Laboratory  was  conducted  in  the  Chukchi  Sea  between  20  July  and 
18  August  from  the  USCGC  STATEN  ISLAND  and  a  drifting  ice  floe.  The  ob¬ 
jective  of  the  study  was  to  obtain  useful  data  on  the  ohenomena  encount¬ 
ered  during  the  previous  summer--namely ,  the  warm  water  intrusion  from 
the  south,  the  resulting  thermal  fine  structure  and  its  effect  on  acous¬ 
tic  transmissions,  and  the  volume  reverberation  produced  by  biological 
scattering  layers.  Reference  1  describes  the  1971  investigations  in  the 
Chukchi  and  Beaufort  Seas. 

A  1-day  helicopter  survey  was  made  in  April  to  observe  spring  ocea¬ 
nographic  conditions.  The  work  in  the  summer  from  the  icebreaker  consisted 
mainly  of  two  cruises,  conducted  during  the  first  and  last  weeks  of  the 
study,  and  was  designed  to  obtain  an  overall  description  of  the  temperature 
and  salinity  distribution  and  the  biological  scattering  layers  in  the  mar¬ 
ginal  ice  zone  (MIZ)  of  the  Chukchi  Sea  as  well  as  some  knowledge  of  the 
changes  that  occur  during  the  summer. 

The  studies  made  from  the  ice  floe  required  a  stable  platform  from 
which  measurements  over  fixed  dimensions  could  be  made.  Temperature- 
depth  profiles  of  the  water  under  the  floe  were  taken  simultaneously 
from  several  stations  to  determine  the  outline  of  the  thermal  layers. 

In  conjunction  with  the  temperature  profiles,  underwater  acoustic  pulses 
were  transmitted  from  one  station  to  another  to  determine  variations  in 
acoustic  propagation  related  to  the  thermal  structure.  The  floe  remained 
in  the  MIZ  until  the  last  day  when  it  was  carried  into  the  Beaufort  Sea 
by  a  strong  current  along  the  coast. 

This  report  contains  the  analysis  of  the  1972  data  which  have  been 
published  in  detail  in  Ref.  2.  Most  of  the  data  required  to  substanti¬ 
ate  the  conclusions  made  in  this  report  are  repeated  here.  The  remaining 
data  in  Ref.  2  may  be  of  use  to  oceanographers  and  arctic  investigators 
with  interests  in  studies  other  than  those  conducted  by  the  authors. 

Local  daylight  time,  Greenwich  time  less  9  hours,  is  used  throughout 
this  report. 
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II.  SUMMARY 


Spring  Conditions 

In  late  April,  a  helicopter  was  employed  to  land  on  various  ice 
floes  from  which  a  temperature  probe  was  lowered  to  the  sea  bottom  to 
obtain  temperature-depth  profiles  of  the  water.  At  the  time,  the  15- 
mile  expanse  of  nearly  open  water  off  Barrow  had  a  temperature  near  the 
freezing  point.  The  only  warmer  water  detected  was  a  -1.4°C  layer  at  a 
depth  of  30  m  near  the  shore-fast  ice  20  miles  southwest  of  Barrow.  The 
origin  of  this  layer  is  unknown. 

Weather 


During  the  August  field  studies  hourly  obseivations  of  the  wind 
speed  and  direction,  atmospheric  pressure,  air  temperature  and  visibil¬ 
ity  were  recorded  manually  while  the  ice  floe  was  occupied.  Floe  ori  -n 
tation  and  water  depth  were  also  noted.  The  weather  ranged  from  clear 
days  with  cold  nights  to  rainy  days  with  relatively  warm  nights  The 
air  temperature  varied  between  32  and  55°F  (0  and  11°C),  the  barometer 
reading  had  a  high  of  30.28  in.  (76.7  cm)  and  a  low  of  29.28  in.  (74  4 
cm),  and  the  wind  speed  ranged  from  0  to  30  mph. 

There  had  been  continual  easterly  winds  during  July  and  it  is 
likely  this  wind  was  the  force  keeping  the  intrusion  of  warm  Pacific 
water  in  the  western  portion  of  the  Chukchi  Sea.  The  sunny  days  and 
light  rains  caused  considerable  surface  melting,  but  this  ablation 
appeared  minor  compared  to  the  loss  of  ice  beneath  the  surface. 


Currents 


e  drift  of  the  occupied  ice  floe,  at  first  slowly  southward  and 
then  rapidly  to  the  northeast,  indicates,  as  in  1971,  a  strong  flow  north¬ 
easterly  along  the  coast  which  appears  to  create  an  eddy  that  carries  the 
ice  lying  north  of  the  stream  southward  into  the  main  flow/current. 


Current  measurements  relative  to  the  floe,  combined  with  floe  drift 
records,  show  that  the  floe  usually  followed  the  main  under-ice  current 
at  60  to  80%  of  the  speed  of  the  water. 


Warm  Intrusions 


In  late  July  the  intrusion  of  Pacific  water  from  the  Bering  Strait 
appeared  as  a  warm  10-m  deep  surface  layer  extending  across  the  western 
portion  of  the  Chukchi  Sea.  This  intrusion  corresponded  to  an  absence 
of  ice  cover  as  far  north  as  latitude  73°. 
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In  the  eastern  Chukchi  Sea  the  intrusion  along  the  coast  of  Alaska 
was  small  compared  to  the  previous  summer  and  was  mainly  on  the  surface. 
The  failure  of  the  warm  intrusion  to  spread  northward  allowed  the  heavy 
ice  cover  to  exist  a  few  miles  north  of  Barrow  in  early  August.  In  mid- 
August  the  intrusion  along  the  coast  reached  8°C  but  remained  within  5 
miles  of  the  coast  at  Barrow. 

Volume  Reverberation 

Extensive  measurements  of  volume  reverberation  were  made  during  two 
1-week  c/uises  in  the  Chukchi  Sea  before  and  after  the  ice  floe  occupa¬ 
tion.  The  instrumentation  consisted  of  two  downward-looking,  narrow- 
beamed,  echo-sounding  systems  at  frequencies  of  38  and  105  kHz.  System 
self-noise  was  such  that  the  instrumentation  could  detect  volume  scatter¬ 
ing  strengths  above  -85  dB  (ref.  1  m-1).  Biological  samples  were  taken 
simultaneously  to  identify  the  scatterers.  In  all,  179  data  files  were 
taken  at  100  different  locations:  98  files  at  105  kHz,  of  which  7  files 
were  taken  by  helicopter,  and  81  files  at  38  kHz. 

Examination  of  the  data  reveals  that  the  acoustic  scattering  strength 
can  vary  greatly  from  location  to  location  even  within  a  few  miles.  Occa¬ 
sional  scattering  layers  were  observed  at  all  depths,  and  many  biological 
species  were  collected.  An  especially  high  scattering  strength  was  ob¬ 
served  in  the  warm  surface  layer  resulting  from  the  intrusion  of  Pacific 
water. 

Scattering  strengths  as  high  as  -40  dB  were  found.  The  scattering 
strength  of  a  typical  layer  was  -50  to  -70  dB.  Numerous  single  targets 
were  also  occasionally  observed  at  all  depths.  The  overall  uncertainty 
in  the  measurement  of  the  scattering  strength  in  a  depth  interval  was 
approximately  ±3  dB.  The  measurements  at  38  kliz  often  correlated  well 
with  those  at  105  kHz,  but  sometimes  they  did  not.  These  data  have  been 
augmented  by  data  collected  in  1973  at  these  frequencies  and  at  60  kHz. 

We  plan  to  improve  our  biological  sampling  methods  and  to  add  a  250-kHz 
system  for  detecting  smaller  organisms. 

Thermal  Fine  Structure 

To  measure  thermal  microstructure,  a  large  ice  floe  was  provided 
with  four  stations  on  the  corners  of  a  140-m  square,  each  station  being 
equipped  for  taking  temperature  profiles.  Profiles  were  taken  simulta¬ 
neously  at  the  four  stations  every  hour,  except  during  the  2  days  when 
profiles  were  taken  at  15-min  intervals  for  6  hours. 

For  the  first  several  days  the  floe  drifted  slowly  southward  through 
the  transition  region  between  the  wi-  "er  water  and  the  warm  intrusion 
along  the  coast.  This  transition  rey. on  contained  thermal  microstructure 
with  temperature  excursions  of  ±0.1°C.  As  the  floe  moved  southward  into 
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the  coastaJ  current,  a  predominant  warm  layer  appeared  at  the  s  irface 
and  larger  microstructure,  ±0.2°  to  ±0.3°C,  appeared  at  the  lower  depths. 
This  structure  represents  horizontal  layers  several  hundred  meters  in 
extent  with  an  average  thickness  of  about  3  m. 

Acoustic  Transmissions 

Transmissions  of  underwater  sound  pulses  at  60  kHz  were  made  between 
two  stations  on  the  ice  floe  at  hourly  intervals  coincident  with  the  tem¬ 
perature  profile  schedule.  Variations  of  5  to  10  dB  were  often  observed 
in  the  region  of  the  thermal  microstructure  at  a  range  of  200  m.  These 
variations  occurred  for  low  angle  transmissions  and  are  predictable  from 
sound  ray  calculations. 

Helicopter  Oceanography 

To  demonstrate  the  feasibility  of  using  helicopters  to  take  measure¬ 
ments  in  the  marginal  ice  zone,  we  often  made  use  of  a  Coast  Guard  heli¬ 
copter.  After  locating  an  open  lead  in  the  ice,  the  helicopter  would 
hover  for  5  to  10  min  just  above  the  water  while  a  transducer  or  temper¬ 
ature  probe  was  lowered  beneath  the  surface.  If  equipment  is  kept  light 
and  measurements  can  be  taken  quickly,  helicopters  provide  a  rapid  means 
of  investigating  a  large  area. 
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II!.  FIELD  TRIP 


The  1-month  summer  operation  consisted  of  two  1-week  cruises  on  ~n 
icebreaker  across  the  Chukchi  Sea  and  a  2-week  drift  on  an  ice  floe.  A 
su.  nary  of  the  activities  is  given  in  Table  I. 
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Itinerary  for  MIZPAC-72. 

Depart  Barrow  on  STATEN  ISLAND 

Barrow  Line 

Wainwright  Line 

Icy  Cape  Line 

Chukchi  Sea  Crossing 

Wainwright  Line 

Search  for  Suitable  Floe 

To  Barrow 

Depart  from  Barrow 

Ice  Camp  Established 

Barrow  Line 

Wainwright  Line 

Wainwright  to  Ice  Camp 

Abandon  Ice  Camp 

Barrow  Line 

Wainwright  Line 

Chukchi  Sea  Crossing 

Southward  Crossing 

Arrive  Nome 


The  USCGC  STATEN  ISLAND  (Figure  1)  supplied  the  main  support  for 
the  operation.  Its  NAVSAT  equipment  provided  valuable  navigational  in¬ 
formation  throughout  the  project  and  ship's  personnel  took  Nansen  casts 
as  requested  to  provide  a  check  on  the  accuracy  of  the  CTD  measurements. 

In  addition  to  providing  reconnaissance  flights,  the  ship's  two 
helicopters  served  as  platforms  for  oceanographic  measurements.  One 
helicopter  contained  a  temperature-depth  probe  with  cable  and  recording 
equipment;  as  the  helicopter  hovered  3  m  above  the  water,  the  probe  was 
lowered  slowly  to  a  depth  of  50  m  and  then  recovered.  The  other  heli¬ 
copter  carried  the  volume  reverberation  equipment;  while  the  helicopter 
hovered,  the  transducer  was  suspended  just  below  the  surface  for  about 
10  minutes  to  take  volume  reverberation  measurements  (see  Figure  2) . 
Before  lowering  any  equipment  into  the  sea,  it  was  necessary  to  drop 
the  helicopter's  sling  cable  into  the  water  to  discharge  the  static 
electricity  that  had  built  up  on  the  helicopter. 
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Figure  1,  STATEN  ISLAND  alongside  ice  floe  camp. 

The  Coast  Guard  radio  station  at  the  Naval  Arctic  Research  Labora¬ 
tory  (NARL)  gave  considerable  assistance  to  the  project,  providing  radio 
communication  between  the  ship  and  personnel  at  NARL. 


Food,  tents,  gasoline,  propane  and  clothing  'equired  for  the  ice  camp 
were  provided  by  NARL.  NARL  also  served  as  a  bas ;  of  operations  and  pro¬ 
vided  food  and  lodging  for  personnel  awaiting  transportation  to  ship  or 
airport.  v 


Instrumentation  aboard  the  icebreaker  for  measuring  temperature  and 
salinity  profiles  and  volume  reverberation  was  provided  by  the  Applied 
Physics  Laboratory.  (This  equipment  has  been  described  in  Ref.  1,  a  re¬ 
port  on  1971  MIZPAC  studies.)  Special  equipment  for  taking  plankton  sam¬ 
ples  was  provided  by  the  University  of  Washington's  Oceanography  Department. 

Two  Omega  receivers  were  prov.'ded  by  the  U.S.  Naval  Oceanographic 
Office.  One  was  operated  at  NARL  '•nd  the  other  was  installed  at  the  ice 
floe  camp.  There  were  various  problems  with  the  Omega  system,  especially 
with  operating  procedures.  For  example,  we  were  not  aware  of  a  1  sec  de¬ 
le.,  between  station  WWV  and  the  start  of  the  Omega  cycle,  nor  that  the 
Japanese  station  on  channel  H  was  operating  on  an  experimental  schedule 
We  were  also  confused  about  the  time  the  New  York  Omega  station  was  re¬ 
placed  by  the  North  Dakota  station  in  segment  D.  In  addition,  sun  flare 
activity  interfered  with  reception  of  both  the  Omega  and  WWV  signals; 
even  when  reception  improved  and  some  of  the  confusing  signals  were  prop¬ 
erly  interpreted,  the  few  fixes  obtained  were  still  a  few  miles  from  the 
NAVSAT  positions. 
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Ice  Floe  Camp 

The  search  for  a  suitable  ice  floe  took  place  on  29  July.  Two  heli¬ 
copters  covered  a  10-mile  square,  stopping  briefly  at  each  corner  of  the 
area  to  obtain  a  temperature  profile.  Several  floes  were  found  and  one 
near  the  ship  was  selected  for  investigation  (see  Figure  3) .  The  command¬ 
ing  officer  of  the  STATEN  ISLAND,  along  with  APL  representatives,  walked 
over  the  floe  and  drilled  several  holes  tc  measure  the  ice  thickness.  The 
large  smooth  area  along  one  side,  an  attractive  feature,  showed  an  ice 
thickness  of  8  ft.  The  ridged  area  was  assumed  to  be  much  thicker.  Because 
the  STATEN  ISLAND  had  to  pick  up  personnel  and  supplies  at  Barrow  before  a 
camp  could  be  established,  we  built  a  10-ft  tower  on  the  highest  point  of 
the  floe,  installed  a  1.35-MHz  radio  beacon,  and  stretched  a  400-ft  antenna 
from  the  tower  to  anchors  in  the  ice. 

On  1  August  the  STATEN  ISLAND  returned  from  Barrow  and  proceeded  to 
the  position  determined  for  the  floe  3  days  earlier.  At  first  the  radio 
direction  finder  gave  some  misleading  information,  but  after  repairs  were 
made  the  floe  was  relocated  easily.  Much  of  the  ice  along  the  outer  edge 
of  the  flat  area  had  vanished  but  the  floe  appeared  to  be  satisfactory  for 
a  2 -week  camp. 

All  building  materials,  equipment  and  supplies  were  transferred  from 
the  ship  to  the  floe.  Four  stations,  at  the  corners  of  a  square,  were 
established  on  the  floe  and  two  APL  representatives  were  assigned  to  each 
station.  Each  station  contained  a  tent  for  living  quarters  and  a  work  hut 
covered  with  a  plastic  tarp  (see  Figure  4).  A  central  tent,  supplemented 
by  a  plastic-covered  hut,  was  used  for  preparing  meals.  Photographs  of 
the  individual  camps  are  shown  in  Figure  5. 

Each  station  was  provided  with  a  temperature  probe,  cable  and  recorder 
for  measuring  temperature  profiles.  Simultaneous  temperature  measurements 
were  scheduled  at  all  four  stations  every  hour.  In  addition,  a  conductiv¬ 
ity  probe  was  included  with  the  profiling  equipment  at  Station  2,  and  a 
transducer  and  receiver  were  provided  at  Station  1  to  record  the  sound 
level  profile  from  a  transmitter  placed  beyond  Station  2.  Equipment  was 
also  provided  for  measuring  wind  speed,  floe  orientation,  water  depth,  air 
pressure  and  temperature,  and  current  magnitude  and  direction. 

Holes  for  the  equipment  were  drilled  with  a  9-in.  ice  auger  or  melted 
with  a  propane- fueled  hot-water  ice  melter  (Ref.  1).  Many  of  the  holes 
passed  through  about  8  ft  of  ice  into  a  2  or  3  ft  pocket  of  water  followed 
by  more  ice.  This  type  of  ice  structure  caused  difficulty  for  the  drill 
and  also  for  the  ice  melter  since  the  hot  water  was  flushed  away  through 
the  cavity. 
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Station  2 
Work  Shelter 


Station  1 
Work  Shelter 


Central  Cook  Tent 


Station  1  Sleeping  Tent 
and  Work  Shelter 


Station  4  Work  Shelter 


Figure  5.  Ice  camp  structures 
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Three  stations  were  in  operation  by  3  August  and  the  fourth  by  5  Aug¬ 
ust.  Profile  data  were  taken  regularly  until  11  August,  when  all  gear  was 
packed  for  loading  on  the  icebreaker  the  next  day. 

Figure  o  is  the  track  of  the  ice  floe  plotted  from  navigational  fixes 
determined  Ky  the  STATEN  ISLAND,  using  NAVSAT  to  determine  the  ship's  posi¬ 
tion  and  then  radar  or  visual  sightings  to  determine  the  relative  position 
of  the  ice  floe.  In  the  figure  a  point  represents  every  6  hours  of  travel 
and  a  short  cross  line  represents  the  start  of  each  day.  The  drift  speed 
increased  to  0.8  kn  on  10  August  and  to  1.3  kn  on  11  August,  causing  consi¬ 
derable  streaming  of  the  profiling  cables.  Fortunately,  most  of  the  planned 
measurements  were  completed  by  this  time. 
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Figure  6.  Drift  track  of  the  occupied  ice  floe. 
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IV.  ENVIRONMENTAL  DATA 


weather  measurements  when  personnel  were  aJaflaMe  CUrrcnt 

Current  Measurements 

T.  Current- measuring  probes  were  lowered  by  hand  from  the  ice  fine 
rhe  measurements  were  time  consuming  and  had  a  low  nrion  w 
data  are  sparse  and  irromiior  a  “  priority,  so  current 

in  Table  II.  “  summary  of  the  measurements  is  given 

Table  II.  Summary  of  current  measurements ,  1972. 

Number  of  Profiles 


Speed  Direction 


3  Aug 

4  Aug 

5  Aug 

6  Aug 

7  Aug 

8  Aug 

9  Aug 

10  Aug 

11  Aug 


5  (20-40  m  only) 


The  equipment  used  to  measure  current  speed  cnncictori  c 

pass,  along  with  a  Vibrotron,  and  was  operated  I  flux-gate  com- 

The  outputs  of  the  compass  and  VibrotroE  were  reL  an^re^Y^1"  line’ 
at  the  surface.  d  and  rec°rded  manually 

ment  stilly  whJeASInrrec^npCtShWe(eHtaken  ^  l0“eri"8  the  i™‘»- 
Uter,  the  instrument  was^held  at  each^epth^Mp"  °UtPut  a"d  r0t0r  count- 
reading.  The  current  measurement  relate  (oT^f,  ”  «el.a  'tatle 
usually  moving.  Direction  profiles  werc^^s”  ts^ent 
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was  stopped  at  every  meter  of  depth  to  record  the  magnetic  direction  of 
the  current  flow.  The  direction  vane  was  then  dragged  along  the  bottom 
to  obtain  the  direction  of  the  floe's  movement.  These  data  are  presented 
in  Figures  7-9.  The  reading  near  the  bottom  is  an  indication  of  the 
actual  speed  of  the  floe  with  respect  to  the  bottom.  When  both  the  speed 
and  direction  of  the  floe's  movement  are  known,  the  current  reading  can 
be  corrected.  6 


Prior  to  7  August  the  current  speed  was  less  than  10  cm/sec  and  the 
direction  was  fairly  constant  at  depths  below  10  m.  On  7  August  the  floe 
picked  up  speed  to  19  cm/sec  in  a  southeasterly  direction  (Figure  6)  and 
the  current  relative  to  the  floe  was  about  15  cm/sec  in  a  northerly  di¬ 
rection  (Figure^9)  .  Adding  the  two  vectors  gives  an  absolute  current  of 
4  cm/sec  at  ,2  true.  Later  in  the  jay  the  relative  current  and  the 
floe  drift  speed  decreased  (Figure  10)  with  an  indication  at  1830  that 
the  current  had  stopped  and  the  floe  was  continuing,  giving  a  relative 
current  to  the  west.  Near  the  end  of  the  day  the  relative  current 
dropped  below  10  cm/sec  in  a  southeasterly  direction.  Since  the  floe 

was  nearly  stationary,  this  indicates  the  speed  and  direction  of  absolute 
current. 


On  8  August  the  floe  was  moving  at  about  8  cm/sec  in  the  same  direc¬ 
tion  as  the  relative  current,  which  dropped  to  5  cm/sec  (see  Figure  11). 


On  9  August  the  floe  picked  up  speed  to  20  cm/sec  at  70° 
a  relative  current  of  15  cm/sec  at  46°  true  (see  Figure  12) 
absolute  current  was  34  cm/sec  at  60°  true.  At  this  time  the 
at  60%  of  the  speed  of  the  water. 


true.  With 
the  computed 
floe  moved 


On  10  August  the  floe  moved  east- lortheast  at  35  cm/sec  and  the  rel¬ 
ative  current  built  up  to  25  cm/sec  ir  the  same  direction,  indicating  an 

absoiute  current  of  60  cm/sec  (see  Figure  13).  Again,  the  floe  moved  at 

oUi  ot  the  water  speed. 

°n  11  August  the  floe  drift  track  shows  a  speed  of  about  100  cm/sec 
A  relative  current  of  20  cm/sec  was  measured  at  15  m,  indicating  an  abso¬ 
lute  current  of  120  cm/sec  if  we  assume  the  direction  was  parallel  to  the 
rift.  At  depths  from  20  to  30  m  the  relative  current  was  near  zero  (see 

Figure  14).  The  increase  in  relative  speed  below  30  m  may  indicate  a 

current  in  the  opposite  direction,  in  which  case  it  would  represent  water 
moving  slower  than  the  floe.  Unfortunately,  no  direction  measurements 
were  made  at  that  time. 
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Figure  7.  Current  measurements  from  the  ice  floe  on  3  and  4  August. 
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Figure  8.  Current  measurements  from  the  ice  floe  on  4  and  5  August, 
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Figure  9. 


Current  measurements  from  the 


ice  floe  on  5-7  August. 
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Figure  10.  Current  measurements  from  the  ice  floe  on  7  August 
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Ftgure  11.  Current  measurements  from  the  ice  floe  on  8  August. 


RELATIVE 
CURRENT 
(CM/SEC) 
SL-O  JO 


I- 


1 


1430 


APL-b'W  7311 


19 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


Figure  12. 
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Current  measurements  from  the  ice  floe  on  9  August. 
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Figure  13.  Current  measurements  from  the  ice  floe 
on  10  August. 


Figure  14.  Current 
measurements  from 
the  ice  floe  on 
11  August. 
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Weather  Measurements 

Weather  observations  and  readings  of  several  instruments  were  taken 
hourly  about  16  hours  each  day  for  the  period  of  the  ice  floe  occupation, 
2-12  August.  Graphs  for  the  10-day  period  are  plotted  in  Figures  15  and 
16. 

The  wind  speed  was  measured  with  a  three-cup  anemometer  mounted  on  a 
12-ft  pole.  A  magnetic  pick-off  sent  electrical  pulses  to  a  counter  at 
the  foot  of  the  pole.  Wind  speeds  were  low  except  for  winds  of  20-30  mph 
on  6  August  and  15-20  mph  on  8  August. 

Wind  direction  was  estimated  oy  facing  the  wind  and  noting  the  di- 
rectior  on  the  compass.  Winds  were  from  the  east  at  first,  changing  to 
westerly  on  5  August  and  continuing  westerly  or  southwesterly  until  the 
floe  was  abandoned  on  12  August. 

The  barometric  pressure  was  recorded  by  an  aircraft  altimeter  with  a 
barometric  correction.  The  correction  cont  -ol  was  varied  to  bring  the 
altitude  to  zero  and  the  barometric  press  was  then  obtained  by  reading 
the  control  setting.  These  data  are  ploti  U  in  Figure  16. 

Air  temperature  was  measured  with  a  thermocouple  meter  sheltered 
from  the  sun.  The  temperature  usually  dropped  to  near  freezing  at  night. 
On  sunny  days  it  rose  as  high  as  52°F  and  during  rainy  days  varied  be¬ 
tween  40°F  and  45°F.  Each  temperature  measurement,  converted  to  degrees 
Celsius,  is  plotted  in  Figure  16. 

The  effect  of  the  wind  on  the  floe  can  be  observed  by  comparing  the 
floe  orientation  and  wind  characteristics  in  Figure  15  with  the  floe 

drift  track  in  Figure  6.  The  high  winds  of  6  and  7  August  from  300°  true 

north  apparently  caused  the  floe  to  move  with  the  wind  toward  120°.  A 
lull  in  the  wind  early  on  8  August  corresponds  with  the  slow  drift  of  the 

floe  on  that  date.  On  9  August  the  floe  apparently  entered  the  coastal 

current  and  was  carried  northeastward  with  steadily  increasing  speed 
despite  the  lull  in  the  wind  for  the  next  two  days. 

This  correlation  between  floe  drift  and  wind  or  current  suggested 
the  possibility  of  predicting  the  drift  track  of  the  floe  from  the  influ¬ 
ence  of  the  wind,  which  was  measured,  and  the  current,  which  had  to  be 
estimated.  The  transfer  of  momentum  from  wind  to  floe  and  from  current 
to  floe  also  had  to  be  estimated.  A  simple  current  pattern  was*  chosen, 
namely,  0.45  kn  at  195°  true  north  except  near  the  coast  where  it  was 
assumed  to  be  1.60  kn  at  070°.  (It  was  assumed  that  the  floe  was  in  the 
first  current  through  7  August,  then  in  the  other  current  for  the  re- 

Sainder  of  the  drift.)  The  assumption  was  made  that  the  drift  velocity 
is  related  to  the  current  C  and  the  wind  $  by  the  vector  equation 
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Figure  16.  Record,  of  pressure  and  air  temperature  at  the  ice  floe  camp 
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5  =  a  £  +  b  W.  The  movement  of  the  floe  was  calculated  for  various  values 
of  a  and  b  and  compared  to  the  measured  movement.  The  best  fit  to  the 
actual  path  was  obtained  with  a  =  1.0  and  b  =  0.1  (see  Figures  6  and  17). 
The  agreement  is  remarkably  good  and  suggests,  but  does  not  prove,  that 
the  assumptions  as  to  the  current  and  the  a  and  b  constants  are  approxi¬ 
mately  correct.  Predictions  based  on  this  calculation  would  fail,  of 
course,  if  a  land  obstruction  were  to  restrict  the  movement  caused  by  the 
wind  and  cause  the  floes  to  pack  together. 


Figure  17.  Predicted  path  of  ice  floe  based  on  wind  and  current. 
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V.  OCEANOGRAPHIC  STUDIES 


The  oceanographic  study  of  the  area  can  be  separated  into  two  parts: 
(1)  the  pattern  of  the  summer  intrusion  of  warm  Pacific  water  through  the 
Bering  Strait  into  the  Chukchi  Sea  and  (2)  the  thermal  fine  structure  pro¬ 
duced  as  this  intrusion  interacts  with  the  existing  cold  and  ice-covered 
water.  The  cruises  during  the  first  and  last  weeks  of  the  project  were 
for  studying  the  intrusion  while  the  main  purpose  of  the  ice  camp  was  to 
study  the  fine  structure. 

Temperature  measurements  were  made  from  a  helicopter  off  the  Alaskan 
coast  in  April  to  determine  the  spring  conditions.  Some  Bering  Strait 
temperature  data  were  obtained  from  the  Arctic  Submarine  Laboratory  for 
comparison  with  the  intrusion  temperatures. 

Most  of  the  temperature  and  salinity  profiling  equipment  was  the  same 
equipment  that  had  been  used  the  previous  summer  (Ref.  1).  A  new  unit  was 
developed,  however,  to  operate  automatically  at  prescribed  time  intervals. 
This  equipment  (see  Figure  18)  lowered  a  temperature-conductivity-depth 
probe  to  the  bottom  at  regular  intervals  and  recorded  the  data  on  paper 
tape.  It  was  used  mainly  for  overnight  operations,  but  should  be  able  to 
perform  for  long  periods  with  improvements  in  tape  storage  and  the  battery 
pack.  A  detailed  description  of  the  automatic  profiler  appears  in  the 
appendix. 
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Spring  Measurements  Off  Barrow  and  at  T-3 

Measurements  were  made  along  the  Barrow  line  on  29  April  1972  to  de¬ 
termine  spring  conditions.  A  helicopter  (Bell  205A-1)  was  rented  from  a 
firm  in  Anchorage  and  provided  with  a  temperature-depth  profiler  and  a 
small  hand  winch  for  taking  Nansen  samples.  The  procedure  was  for  the 
helicopter  to  land  on  the  ice  and  remain  with  rotor  turning  while  the 
probe  was  carried  to  the  ice  edge  and  lowered  to  the  bottom  or  to  the 
full  cable  length  (60  meters) .  The  recording  instrumentation  remained 
aboard  the  helicopter. 

The  first  stop  was  made  at  the  edge  of  the  shore-fast  ice  about  a 
mile  from  shore.  Stops  were  planned  at  10-mile  intervals  but  there  was 
so  much  open  water  and  thin  ice  that  the  second  station  was  at  20  miles. 
Succeeding  stations  were  spaced  every  10  miles  out  to  70  miles. 

A  second  series  of  four  measurements  at  20-mile  intervals  was  taken 
at  the  edge  of  the  shore-fast  ice  along  the  coast  toward  Wainwright.  The 
last  station  of  this  series  was  midway  between  Wainwright  and  Icy  Cape. 

Nansen  samples  and  reversing  thermometer  readings  were  taken  at 
three  of  the  stations,  two  in  the  first  series  and  one  in  the  second. 

One  thermometer  failed  but  the  other  showed  fair  agreement  with  the 
thermistor  used  for  profiling. 

The  temperature  profiles  are  shown  in  Figures  19  and  20.  In  these 
profiles  the  temperature  scale  is  greatly  expanded  to  bring  out  the  fine 
structure  in  the  nearly  isothermal  water.  Reversing  thermometer  readings 
are  shown  as  dots  next  to  the  profile  and  salinities  are  written  near  the 
dots. 


Salinities  were  about  33%o  60  miles  off  shore.  The  high  salinity 
value  at  5  m  in  the  first  profile  would  indicate  a  density  inversion 
and  may  be  in  error. 

The  water  temperature  throughout  the  area  was  only  slightly  above 
the  freezing  point  for  water  of  this  srlinity.  Some  layering  appears  in 
the  profiles  taken  at  Barrow  and  20  mil'js  down  the  coast,  but  the  measure¬ 
ments  show  that  there  was  no  appreciable  amount  of  warm  water  entering  the 
area  at  that  time.  The  considerable  amount  of  open  water  off  Barrow  was 
not  related  to  the  presence  of  a  warm  intrusion,  since  layers  only  0.1°C 
warmer  than  the  main  body  of  water  hardly  could  be  considered  significant. 

A  series  of  measurements  at  Fletcher's  Ice  Island  (T-3)  was  taken  in 
May,  One  of  these  has  been  plotted  and  included  here  for  comparison  (see 
Figure  21).  The  temperature  remained  near  the  freezing  point  to  a  depth 
of  50  m,  below  which  both  temperature  and  salinity  increased. 


Figure 


Figure  20. 
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Oceanographic  measurements  along  the  coast  southwest  of  Barrow. 
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Figure  21.  Profiles  at  T-3  in  May  1972. 


Chukchi  Sea  Intrusion 

Two  one-week  cruises  were  planned  to  investigate  temperature  and 

Slnii£  c™dltl0ns  nlon8  the  marginal  ice  zone  (MIZ)  across  the  Chukchi 
bea.  The  first  started  at  Barrow  on  20  July  (see  Figure  22)  and  included 
lines  across  the  MIZ  at  Barrow,  Wainwright  and  Icy  Cape.  After  a  west¬ 
ward  cruise  nearly  to  Herald  Island  and  back,  a  final  crossing  of  the  MIZ 
was  made  at  Wainwright.  The  second  cruise  started  on  13  August  and 
included  lines  at  Barrow  and  Wainwright.  The  westward  cruise  was  then 
made  further  north  than  before  in  order  to  remain  close  to  the  receding 

tCA  !jge',  a Z  ab°Ut  73-N  169° W  the  temPerature  measuring  equipment  failed, 
and  the  study  was  terminated.  Between  these  cruises  a  few  short  lines 

were  run  to  observe  conditions  in  the  vicinity  of  Wainwright  and  Barrow. 
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The  temperature  and  salinity  profiles  for  each  of  the  stations  shown 
in  Figure  22  are  presented  in  Ref.  2.  These  data  have  been  used  to  con¬ 
struct  sectional  views  of  the  temperature  and  salinity  conditions  at  each 
line  of  stations  across  the  MIZ. 

Barrow  Line,  20-21  July  1972 

Representative  profiles  along  the  Barrow  Line  are  presented  in  Figure 
23.  Sectional  views  of  the  temperature  and  salinity  values  along  this 
line  are  presented  in  Figure  24.  A  slight  warming  and  a  surface  layer  of 
low  salinity  water  can  be  noted  near  the  coast.  Two  days  later  a  series 
of  stations  from  a  helicopter  gave  similar  results  (see  Figure  25). 

Wainwright  Line,  24  July  1972 

Sample  profiles,  along  the  Wainwright  line  are  shown  in  Figure  26. 
Sectional  views  of  temperature  and  salinity  are  shown  in  Figure  27.  The 
intrusion  is  10  m  thick  with  a  maximum  temperature  of  5°C  and  a  salinity 
of  about  29%0. 

Icy  Cape  Line,  25  July  1972 

The  large  surface  layer  shown  for  this  region  in  Figures  28  and  29 
increases  in  depth  and  temperature  as  one  proceeds  southwest  toward  the 
Bering  Strait.  The  maximum  temperature  is  8°C. 

Westward  Crossing,  25-27  July  1972 

Representative  temperature  and  salinity  profiles  for  this  crossing 
of  the  Chukchi  Sea  are  shown  in  Figure  30.  Sectional  views  are  presented 
in  Figure  31.  It  appears  that  the  intrusion  has  split  into  two  parts 
separated  by  an  ice-covered  region.  The  salinity  of  the  intrusion  is  not 
much  lower  than  that  of  the  existing  water.  A  series  of  temperature  pro¬ 
files  from  a  helicopter  on  a  line  to  the  north  (Figure  32)  shows  a  tongue 
of  4°C  water  extending  into  the  cold  region.  No  salinity  data  were  taken 
along  this  line  because  the  lightweight  profiler  used  in  the  helicopter 
did  not  include  a  conductivity  cell. 

Wainwright  Line,  28  July  1972 

At  the  end  of  the  first  cruise  the  Wainwright  line  was  repeated  to 
determine  what  changes  were  taking  place.  The  sectional  temperature  dia¬ 
gram  of  Figure  33  shows  that  the  intrusion  is  similar  in  size  to  that  on 
24  July  (Figure  27)  but  is  now  well  away  from  the  coast.  A  continuation 
of  the  measurements  to  the  northeast  (Figure  34)  shows  a  pool  of  consid¬ 
erably  warmer  water  which  is  somewhat  cut  off  from  the  rest  of  the  intru¬ 
sion  by  colder  water  to  the  southwest.  The  warmer  water  may,  however, 
connect  to  the  main  body  of  the  intrusion  directly  west  of  this  pool. 
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Figure  23.  Temperature  profiles  along  the  Barrow  line. 
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Figure  24.  Isotherms  and  iso¬ 
halines  for  a  section  along 
the  Barrow  line. 
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Figure  25.  Isotherms  for  a  section  along  the  Barrow  line. 


*i 


TEMP  (*C) 


TEMPrci 


24  JUL  72 


Figure  26.  Temperature  profiles  along  the  Wainuright  line. 
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Figure  27,  Isotherms  and  Isohalines  for  a  section  along  the  Wainwright  line. 


36 


APL-UW  7311 


DEPTH (METEPS) 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


TEMPfC) 


30  32 

SALINITY  (•/..) 


40r 

165*  W  LONG 
2010  25  JUL 


5Qt- 


TEMP  (*C| 


|T«*W  LONG 
0155  IT  JUL 


FIRST  CROSSING 
CHUKCHI  SEA 


Figure  30.  Temperature  profiles  for  the  cruise  across  the  Chukchi  Sea 
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Figure  31.  Isotherms  and  isohalines  for  an  east-west  section  aoross  the 
Chukchi  Sea. 
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Figure  33.  Isotherms  for  a  section  along  the  Wainuright  line. 
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Figure  34.  Isotherms  for  a  section  along  a  line  northeasterly  from  the 
end  of  the  Wainwright  line. 


Traverses,  4-7  August  1972 

During  the  ice  camp  occupancy,  the  icebreaker  was  available  for 
three  traverses  of  the  area  between  Wainwright  and  Barrow.  The  4  August 
traverse  in  Figure  35  again  shows  a  split  in  the  intrusion,  one  part 
along  the  coast  and  one  at  50  miles  out.  The  5  August  traverse,  Figure 
36,  which  ends  at  about  the  center  of  the  second  Wainwright  line  (Figure 
34),  gives  the  impression  that  the  intrusion  has  been  forced  to  the  coast. 
The  third  traverse  from  Wainwright  northward  then  eastward  toward  Barrow, 
Figure  37,  shows  a  concentration  of  the  intrusion  along  the  coast  at 
Wainwright  and  very  little  22  miles  off  Harrow. 

Barrow  Line,  13  August  1972 

The  second  one-week  cruise  started  after  the  ice  camp  was  abandoned. 
The  profiles  of  Figure  38  and  the  sections  presented  in  Figure  39  show 
that  the  intrusion  is  still  crowding  the  coast,  with  only  remnants  beyond 
8  miles. 


APL-UW  7311 


41 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


OtSTANCE  fROM  BARROW  (MILES) 


Figure  35.  Isotherms  for  a  section  along  a  line  west  from  Barrow. 
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Figure  36.  Isotherms  for  a  section  along  a  line  north  from  Wainwright. 
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Hgvre  38.  Temperature  profiles  along  the  Barrow  line. 
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Figure  39.  Isotherms  and  isohalines  for  a  section  along  the  Barrow  line. 
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Wainwright  Li  le,  14  August  1972 

The  highest  temperature  of  the  intrusion,  over  8°C,  is  shown  in  the 
profiies  taken  near  Wainwright  on  14  August  (see  Figure  40).  The  extent 
of  the  warm  high-salini  y  surface  layer  can  be  seen  in  the  sections  in 

ii8th6  i1'  lt  appears  t,  at  the  surface  water  has  moved  toward  the  coast 
in  the  same  manner  as  the  manned  ice  floe. 
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Figure  40.  Temperature  profiles  along  the  Wainwright  line. 


Second  Chukchi  Crossing,  15-16  August  1972 


,se^s  is"  "8“hr: 

&5T"  t0  be  a"  iCe  pr0je,:ti0n  fr°"  ‘he  "°«h  that  resists  the'"8’ 

three^rfiiir^tx:^  ss* 

sr*s:  sd  s 

la^^rj^iirc^^StmerS'a  In  the  SeCOnd  fiEUre 

coider  water  Era.  the  “4 


46 


APL-UW  7311 


APL-UIV  7311 


4 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


to*  ret 


n  lat 

i 


•  KM 

A  lAf 


SfCOWO  CftOSSMG 
CMUHCHI  SCA 
IS-*  AUG  72 


Figure  42. 


Tempevatw*e  profiles  for  the  second  cruise  across  the 
Chukchi  Sea. 


5ssr»!?^sri5 “■»»-*  - 

The  Arctic  Submarine  Laboratory  provided  some  Berine  Strait  surface 
layer  temperatures  taken  in  shallow  water  near  Wales  AUska  ke  jfn 

-  =r. siSsre? 

SK’4*".  SSs 

station  at  Wales  is  out  of  the  main  stream  and  the  warm  intrusion  through 

off  wllef  Sreral  dayS  earlier  than  the  observed  temperature  increase 
wf!es.  The  temperature  peaks  at  Wales  were  1  or  2°C  lower  in  1971 

thecas tern^huk chi  SeTin  ^  tCmP“  °f  the  ^ -  in 

e-fr 

example,  an  increase  in  temperature  of  8°C  will  cause  sound  rays  less 

If  thl\  6  horizontal  to  be  bent  downward.  Examples  of  the  effect 
the  layer  on  sound  transmission  are  given  in  Figure  48  Tn  tho 
sound  ray  plot  the  source  is  in  the  theLcuJe;  in  the  LeJ  p^t 
source  io  below  the  thermocline  and  the  9°  ray  is  refracted  doSnward^ 
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Figure  45.  Isotherms  cf  maximum  temperature  for  the  second  cruise  across  the  Chukchi 
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Figure  46.  Isotherms  of  maximum  temperature  for  August  1971. 
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Bering  Strait  temperature  as  recorded  near  Males  by  the 
Arctic  Submarine  Laboratory. 
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Thermal  Microstructure 


Measurements  taken  simultaneously  at  two  stations  the  year  befo 
had  shown  that  the  thermal  microstructure  was  distributed  m  layers  of 
considerable  horizontal  extent.  In  order  to  obtain  a  better  understand¬ 
ing  of  Se  size  and  movement  of  the  layers,  four  stations  were  established 
on8 the  floe  at  the  corners  of  a  140-m  square  and  provided  ":th 
for  recording  temperature  profiles  at  regular  intervals.  Carefully  call- 
bratec^thermistors  ensured  that  the  profiles  at  the  four  statrons  could  be 

compared . 

The  profile  data  were  transferred  from  paper  tape  to  magnetic  tape 
for  processing  by  the  CDC  6400  computer  at  the  University  of  Washington. 

A  special  program  was  developed  to  read  the  tapes,  calculate  and  p*°*  the 
nrofiles  select  and  characterize  all  thermal  layers,  and  punch  a  card 
for  each* layer.  The  profiles  were  then  examined  and  compared  to  determine 
Zch layers  were  identifiable  in  profiles  taken  at  other  times  or  at 
other  stations.  Such  layers  were  numbered  and  all  cards  for  that  layer 

were  segregated. 

To  facilitate  computations,  a  laver  was  defined  as  what  might  better 
be  described  as  a  "half  layer"  (Figure  49).  For  many  closely  fPaoed 
lavers  the  definition  is  a  convenient  one.  A  negative  layer  is  shown 
i/the’ figure;  the  next  layer  below  it^would  be  positive.  Layers  with 
temperature  differences  less  than  0.03  C  were  ignored. 
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Figure  49.  Authors'  definition  of  a  thermal  layer. 
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^h®  first  measurements  of  the  temperature  nrnfi  io  a*  i  •  j- 

a  fairly  smooth  transition  from  0«C  a^™^  ^^75^  lo  T* 
However,  when  the  profiles  were  replotted  at  an  emanHeH  croi 
ture  with  variations  up  to  0.2«C  was  ob^rvod  ?o  aTplh 

occupation^'  ^  Character  hourl>'  but  Parted  throughout  the  ice  fl« 


The  most  complete  series  of  hourly  temperature  profiles  was  made  at 
Camp  1.  Figure  50  is  a  plot  of  this  seriesfor  the  8-day  period  l" 
the  figure  each  successive  profile  is  offset  0.5°  from  tiTprevious  one. 

From  4  to  9  August  an  attempt  was  made  to  take  hourly  profiles  at 

u  ianyroSt;hreeSstaJtan60USly-  6eCaUSe  °f  e<tuiP>e” ‘  failure,  howler, 
for  comparisoJ  iTFiguie  5l.Wer6  °Peratin8'  Th««  P™*^  are  plotted 

On  8  and  9  August  profiles  were  taken  simultaneously  at  all  fn„r  ct* 

o^i  ttzrstt?  r  Hs  Hr 

temperature  “«££.*(&  Hnuy""  IZ^VrV/e  ^  "T 
on  each  day  was  calculated  a ■  ?Dsefved  for  the  6-hour  period 

anH  tw  y  calculated.  The  distribution  in  temperature  difference 

s.’ »“  ““  ariz:1:. 


Sharpness  (degrees/meter) 
Positive  Layers  Negative  Layers 


8  Aug 

9  Aug 


0.092 

0.080 


0.057 

0.107 


The  thermal  layers  can  usually  be  identified  by  their  temperature  denth 
and  shape  on  the  temperature-depth  profile  Manv  nf  V  ,pth 

£  x:  se  3^  fc  i :H 
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Figure  50.  Hourly  temperature  profiles  at  Camp  1 
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Figure  50.  Hourly  temperature  profiles  at  Camp  1 ,  oont. 
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Figure  50.  Hourly  temperature  profiles  at  Camp  .\3  oont 
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Fxgure  50.  Hourly  temperature  profiles  at  Camp  1,  oont. 
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Figure  51.  Simultaneous  temperature  profiles  at  four  stations,  oont 
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Figure  51.  Simultaneous  temperature  profiles  at  four  stations t  oont. 
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Figure  51.  Simultaneous  temperature  profiles  at  four  stations }  oont. 


APL-UW  7311  69 


DEPTH  (METERS)  DEPTH  (METERS)  DEPTH  (METERS) 


i 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 
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Figure  SI.  Simultaneous  temperature  profiles  at  four  stations. 
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Figure  51.  Simultaneous  temperature  profiles  at  four  stations,  oont. 
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Figure  52.  Temperature  profiles  for  Camps  1-4 ,  8  August  1972 
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Fi9Ure  5 4'  Di8trttution  in  temperature  difference  of  thermal  layers. 


*ppear  *?  Very  irregular  and  attempts  to  determine  the 
h  pe  and  size  of  an  individual  layer  have  been  unsuccessful  Therefr.™ 
we  havf-  concentrated  on  studying  the  layers  *  ’j  Therefore 

~  2aracterlstics  •  ^ 

-J'srsL.iiS:  for  — »ie- 8  “■> » Um. 

On  8  August  the  floe  moved  very  slowly  0  15  kn  at  Tin0  u 
The  relative  current  was  0.14  kn  at  146°  t^ue  with  an  8 

tChe\nLfat  ;he,"atf.Carried  the  ™  ™o"thea“erl'Jnd5"ctIo„  bui  1_ 

the  head  wind  slewed  its  progress.  On  the  average,  each  layer  ras  id^r 

tive^curr  T*  'UCCessive  Profiles  (135  minutes)  Considering  the  rela-’ 
tH  Jt  L  assuming  that  the  floe  was  moving  past  layers  fixed  in 
the  water,  the  average  layer  length  would  then  be  600  m. 


On  9  August,  the  floe  moved  much  faster,  0.42  kn  at 


96°  true,  with 


a  relative  current  of  0.36  kn  at  41°  tm?  a  t  „  • ,  .  ,  ,  - 

little  effort  r>n  *  .  '  A  3_kn  tai1  Wlnd  probably  had 

r‘utiva2™!Sof  1300"“"®  Pr°' 
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R/D 


Figure  55.  Number  of  stations  that  Would  observe  random  layers  with  a 
ratio  of  layer  radius  to  station  spacing  of  R/D. 


These  calculations  assume  that  the  floe  is  moving  over  fixed  pre¬ 
existing  layers  The  large  difference  in  layer  size  (600  m  vs  1300  m) 

assumPt^on  15>  incorrect.  On  the  other  hand,  a  hypo- 
esis  that  the  layers  are  formed  m  the  vicinity,  decay  slowly  because 
of  mixing,  and  should  therefore  be  observed  for  a  fairly  constant  time 
(assuming  the  floe  stayed  in  the  vicinity  of  the  layer  long  enough)  would 
mean  that  the  layers,  which  were  observed  for  about  two  hours  on^oth 
days,  were  larger  than  1300  m.  The  true  picture  is  perhaps  somewhere  be¬ 
tween  these  two  hypotheses--some  layers  deteriorate  and  others  move  out 
or  view. 


Another  indication  of  layer  extent  can  be  obtained  by  assuming  that 
the  temperature  difference  of  a  layer  de-reases  linearly  with  distance 
By  comparing  the  observed  difference  at  two  camps  one  can  compute  Se  dis- 
tance  to  the  end  of  the  layer.  Doubling  this  distance  gives  a  rough  esS- 
mate  of  the  horizontal  extent  of  the  layer.  This  procedure  yields  layeJ 
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lengths  starting  at  1000  m  and  scattering  to  three  times  that  size  for 
both  dates. 

Considering  the  spread  in  the  various  estimates  of  layer  size  made 
above,  it  appears  that  the  observed  layers  were  very  irregular  in  shape 
and  varied  widely  in  size,  from  about  250  m  to  1500  m.  In  this  study  no 
noticeable  correlation  between  size  and  any  other  characteristic  of  the 
layer  or  of  the  MIZ  was  detected. 
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VI.  ACOUSTIC  PROPAGATION 


The  short  range  transmission  of  acoustic  pulses  in  the  10-100  kHz 
range  is  important  for  submarine  operation  in  the  marginal  ice  zone.  The 
prevalent  thermal  microstructure,  especially  when  it  consists  of  horizontal 
layering,  results  in  considerable  distortion  to  low  angle  transmissions. 

In  order  to  measure  and  study  the  effect  of  the  thermal  layers,  a  transmit¬ 
ter  was  installed  at  one  station  on  the  ice  floe  and  a  receiver  at  another. 
The  transmitting  transducer  was  held  at  a  fixed  depth  while  the  receiving 
transducer  was  raised  or  lowered  to  obtain  a  sound  level -depth  profile. 

At  the  same  time  a  thermistor  and  Vibrotron  pressure  transducer  attached 
just  above  the  receiving  transducer  obtained  a  temperature  profile.  These 
measurements  were  taken  hourly  from  4  to  9  August.  During  the  last  two 
days  on  the  floe  some  additional  measurements  were  taken  at  other  depths 
and  ranges . 

Instrumentation 

The  acoustic  transmission  equipment  consisted  of  an  electronic  trans¬ 
mitter,  aa  electronic  receiver,  two  identical  ring-stacked  transducers  (each 
on  60  m  of  cable),  and  three  precision  pulsers.  The  receiver  output  had 
a  frequency  proportional  to  the  Log  of  the  input  level  and  was  fed  to  a 
paper-tape  punch.  An  oscilloscope  was  used  to  monitor  the  receiver.  A 
block  diagram  of  the  acoustic  transmission  equipment  is  shown  in  Figure  56. 

The  precision  pulsers  were  portable  units  which  contained  their  own 
batteries,  used  oven-stabilized  crystal  oscillators,  and  had  rated  drift 
rates  of  less  than  one  part  in  10®  per  hour  after  a  1-hour  warmup.  In 
practice,  the  difference  between  two  units  was  less  than  200  ysec  in  a 
24-hour  period.  In  order  to  obtain  2-day  operation  of  the  pulsers,  ex¬ 
ternal  lead  acid  batteries  were  used  at  the  receiving  and  transmitting 
units.  The  third  pulser  was  used  to  carry  the  synchronization  from  one 
pulser  to  the  other  at  the  initiation  of  the  experiment  and  when  batteries 
were  changed. 

The  transmitter's  pulse  width  was  4  msec  and  the  output  voltage  was 
44  V  rms ,  giving  an  acoustic  output  of  85  dB  per  volt  (ref.  1  ybar  at  1 
m)  in  the  center  of  the  transducer  beam. 

A  block  diagram  of  the  receiver  is  shown  in  Figure  57.  The  dynamic 
range  of  the  input  voltage  was  -100  dBV  to  -60  dBV  (i.e.,  10  yV  to  1  mV) 
with  the  attenuator  set  at  zero.  The  attenuator  range  was  38  dB  in  2-dB 
steps,  and  thus  the  full  range  of  the  receiver  was  -100  dBV  to  -22  dBV. 

The  0  dB  setting  was  used  exclusively  for  all  profiles  taken  on  this  trip. 
The  output  of  the  amplifier  was  sampled  for  1  msec  at  a  time  determined 
by  a  delayed  gate.  The  peak  value  was  held  and  applied  to  the  logarithmic 
amplifier.  The  output  of  this  amplifier,  presented  on  a  panel  meter  for 
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Figure  56.  Block  diagram  of  the  acoustic  transmission  equipment 


operated  i/the  range^f  12-25^2°  aTJ'“lta8e‘'arjabl«  oscillator  which 
summed  with  the  gating  pulse  to  provide  an^,?!  the  a”pIifi  er  was  also 
that  the  gate  could  4  L,«,tS 
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The  transducers  were  60-kHz,  ring-stacked  magnetostrictive  units  hav¬ 
ing  omnidirectional  (within  ±0.5  dB)  patterns  in  the  horizontal  plane  and 
a  26  degree  pattern  (half-power  points)  in  the  vertical  plane.  Their  open 
circuit  receiving  sensitivity  was  -99  dBV/ybar;  their  transmitting  output 
was  52  dB  per  volt  (ref.  1  ybar  at  1  m) ;  and  their  impedance  at  the  cable 
connector  was  80  +  jO  ohms.  The  transducers  were  identical  within  a  small 
tolerance  and  were  used  interchangeably  for  transmitting  and  receiving. 

Data  recording  instrumentation  included  a  data  signal  converter  and 
a  paper-punch  unit  (see  Figure  58) .  The  signal  converters  used  for  the 
TD  and  CTD  systems  were  capable  of  converting  each  of  three  frequency- 
modulated  input  signals  to  four-character  BCD  numbers.  The  converters 
utilized  a  BCD  counter  that  stores  the  number  of  1-MHz  clock  pulses  gene¬ 
rated  during  255  complete  cycles  of  the  sensor  processor  output  signal. 

The  output  of  the  BCD  counter  is  gated  by  the  functional  timing  controls 
and  at  the  appropriate  time  is  supplied  to  the  paper  tape  unit  for  record¬ 
ing.  The  input  sensor  signals  are  sequentially  multiplexed  to  produce  a 
data  sample  that  consists  of  three  BCD  outputs,  a  time  code-word  used  to 
represent  the  relative  time  of  the  sample  from  the  start  of  the  drop,  and 
an  externally  set  12-bit  identification  code.  A  data  sample  was  obtained 
in  approximately  0.8  sec 

The  punched-paper-tape  output  for  each  profile  was  about  25  ft  long. 
The  tape  was  usually  fed  into  a  box  and  rolled  up  after  lowering  the 
probe.  Although  it  was  necessary  to  keep  the  tape  dry,  it  was  possible 
to  operate  the  punch  in  the  open  on  many  days. 


Figure  58.  Block  diagram  of  the  data  recording  instrumentation. 
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nun r®ductl0n  system  was  provided  in  the  field  to  convert  the 
p  nched-tape  data  to  a  plot  of  temperature  and  salinity  versus  depth 
OccasionaHy  sound  velocity  and  density  were  also  included.  This  sys- 
tem  consisted  of  a  Hewlett-Packard  tape  reader,  desk  calculator,  printer 
and  plotter  A  50-m  profile  was  plotted  in  about  15  minutes.  The  calcu- 
lator  rapidly  performed  the  lengthy  calculations  of  salinity,  sound  speed 
and  density.  The  system  was  quite  flexible  since  the  plotter  could  be 
set  to  plot  each  variable  at  any  desired  scale.  The  plotter  was  the 
slowest  link  m  the  system,  with  the  tape  reader  only  slightly  faster. 


Transmission  Measurements 

Underwater  transmissions  were  made  from  one  station  on  the  ice  floe 
to  an°  her  to  study  the  effect  of  the  ther[nal  microstructure“ne  J“rtflM 

range  transmissions.  The  transmitter  was  set  at  a  depth  of  33  m  about 
00  m  away  from  Camp  1  and  on  a  line  extending  beyond  Camp  2  so  that  the 

^ofilerSitnrPath9C°llad  1°  mCnitored  by  both  the  temperature-salinity 
Camp  2  and  the  temPerature  profiler  at  Camp  1.  A  receiving 
l,vPidUCeriaS  attached  to  the  temperature  probe  at  Camp  1  and  sound 
level  profi.es  were  measured  routinely  at  the  camp  at  the  same  time  as 

d^rinrtheTdavrthe165'-  ^ l  2°°  intensity  Profiles  were  measured 
auring  the  8  days  the  equipment  was  installed  on  the  ice  floe  Each  n™ 

i  e  consisted  of  sound  level  measurements  at  200  successive  depths  as 
^rreeC^i;tngtLrtan:nU2edr  °P  A  pr°file  taken lach 

of  15  minutes  /S  3  Series  W3S  °btained  with  an  interval 


An  example  of  a  large  variation  in  the  received  acoustic  intend  tv 

illKTS: WblCh  Sh°WS  the  S3linlty  and  te-pemi™5^. 

es  at  each  end  of  the  transmission  path,  and  the  intensity  nrofile  at 
the  receiver  In  this  case  the  large  anomaly  in  the  Sens ^profile 
can  be  related  to  a  warm  layer  which  bends  the  sound  down  causing  1 

For^n  r  Jn^nsity  at  the  dePth  of  the  layer  and  an  increase  below  it 
For  most  of  the  measurements  the  thermal  layers  were  smaller  and  mnrl 

Jigure°U60)?nd  ^  intensitr  Pr°files  -re  more  complex  (see 


The  transmission  measurements  taken  during  the  8-day  period  are 

series  taken  onTlndYi®  ^  rep°rt  (Ref>  2)  *  two  6-hour 

series  taken  on  8  and  9  August  are  of  special  interest  because  of  the 

short  interval  (15  min)  between  measurements.  These  sound  levpl  nrn''-i 

are  plotted  in  figures  61  and  62  alongside  the  cSrrLpZing  a^iL  5  $ 

sound  velocity  profiles  at  Ca„p  2.  The  decrease  in  sound  ?evet  Hucw 

.tions  curing  the  afternoons  of  both  days  is  easily  attributed  to  the 

decrease  m  the  amount  of  thermal  structure  at  that  time. 
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Figure  60. 


Sound  intensity  variation  resulting  from  several  small 
thermal  layers. 
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Figure  62b.  Intensity  profiles  at  Camp  1,  9  August  1972. 
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Effect  of  Layers 

If  the  layers  extend  throughout  the  transmission  path  with  little 
change  in  shape,  a  set  of  sound  rays  can  be  readily  computed  along  with 
the  intensity  levels  along  the  ray  paths.  A  ray  trace  computer  program 
obtained  from  Naval  Undersea  Center,  San  Diego,  (Ref.  3)  was  used  to  cal- 

tHe  ray  PathS  and  tHe  S°Und  leVel  al0ng  each  raX-  An  additional 
forPpn^  pr0gram  was  written  t0  determine  by  interpolation  the  sound  level 
f°  ®  Ch  ray  at  a  Prescribed  range  and  to  plot  sound  level  versus  depth 
at  this  range.  These  calculated  profiles  are  compared  to  the  measured 

with1  the  h"  fFlgUle  ?3‘-  TH!  measurements  selected  for  comparison  were  those 
ah*h,th  beft  ^relation  between  the  temperature  profile  at  Station  2 

profiled  Station  ^he,'/ay  along  the  transmission  path,  and  the  temperature 
t  Station  .,  ,t  the  receiving  end.  Profiles  with  pronounced  layers 
near  the  transmitter  jepth  were  also  favored. 

n™n!fe/t\SUltSi.infCate  that  the  intensity  variations  can  be  roughly 
predicted  through  the  use  of  sound  ray  computations.  The  failure  of  the 
computations  to  predict  the  variations  in  detail  is  assumed  to  result  from 
a  change  in  the  sound  velocity  profiles  along  the  transmission  path. 

Examination  of  the  measured  sound  level  profiles  shows  that  there  is 

frin\hn  unpredl^ted  3  t0  5  dB  fluctuation  at  depths  considerably  removed 
from  the  source  depth.  These  fluctuations  could  be  caused  by  a  wobble  of 
-  in  the  transducer,  but  the  short-range  measurements  show  a  variation 

fipuration°UtTh  f  ^  S?**  *e"ion) »  indicatinS  a  stable  transducer  con- 
figuratron  These  unpredicted  fluctuations  may  be  caused  by  horizontal 

bending  of  the  sound  rays  and  the  resulting  interference. 

S°J"e  of  the  measured  sound  level  profiles  show  violent  fluctuations 
at  depths  near  the  transmitter  depth.  These  fluctuations  are  usually  in 
a  region  where  two  or  more  ray  paths  intersect.  If  rays  intersect  at 
ngle  P  (see  Figure  64)  and  are  in  phase  at  A,  they  will  also  be  in  phase 
we  havf  S  15  3  mUltiple  °f  a  ^velength.  For  s  equal  to  one  wavelength, 


sin  P  =  d. 

h 


from  which 


h  =  s  esc  P 
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VELOCITY  PROFILE  1M/SEC) 


SOUNO  LEVEL  (DB) 


!>  measured 

\  - 


8J  SJ 

5  AUG  1972  AT  1800  RANGE  =  239  YD 


VELOCITY  PROFILE  1  M/SEC  1 

1438  1440  1442  1444 


SOUNO  LEVEL  (DB) 

40  SO  60 


;  PREDICTED 


ry 


6  AUG  1972  AT  1400 


;  ~v 

T 


RANGE  =  239  YD 


Figure  63.  Comparison  of  measured  to  predicted  sound  evel  profiler. 
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Figure  63. 


Comparison  of  measured  to  predicted  sound  level  profiles ,  cont 
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7  AUG  1972  AT  2300  RANGE  =  219  YD 


SJ 

3  AUG  1972  AT  0100  RANGE  =  220  YD 


Figure  62.  Comparison  of  measured  to  predicted  sound  level  profiles t  cont 
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VELOCITY  PROFILE  (M/SEC) 


SOUND  LEVEL  (dB) 


•0 

RANGE  =  224  YD 


Figure  63.  Comparison  of  measured  to  predicted  sound  level  profiles 3  cont 
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Values  of  h  for  small  angles  of  intersection  are  given  in  Figure  64.  In 
regions  where  rays  intersect,  the  sound  level  should  fluctuate  through  a 
cycle  when  the  receiver  is  moved  vertically  from  A  to  B.  If  the  two  rays 
are  of  equal  amplitude  the  receiver  would  record  a  sound  level  with  a 
peak  of  +6  dB  when  in  phase  and  a  minimum  value  representing  noise  when 
180°  out  of  phase. 


I 


Figure  64.  Intensity  fluctuation  when  rays  intersect 
(wave  front  curvature  exaggerated) . 
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Some  of  the  sound  ray  diagrams  were  checked  for  intersecting  rays. 
Where  rays  intersected  there  was  usually  a  fluctuation  with  depth  witha 
period  somewhat  similar  to  that  given  by  the  equation  for  h  above. 

Sound  level  measurements  were  made  for  several  different  source 
depths  and  ranges  to  obtain  some  indication  of  the  variation  of  intensity 
fluctuations  with  range.  These  measurements  are  plotted  in  Figures  65, 

66,  and  67  along  with  the  sound  level  profile  predicted  from  sound  ray’ 
diagrams  for  the  measured  sound  velocity  profile.  The  sound  velocity 
profile  varied  so  much  during  the  5  hours  required  to  complete  the  mea¬ 
surements  that  an  accurate  range  dependence  could  not  be  determined.  (A 
very  small  change  in  the  sound  velocity  profile  near  the  transmitter  will 
produce  large  changes  in  the  sound  level  profile  several  hundred  meters 
away.)  However,  it  appears  that  the  very  large  fluctuations  that  large, 
sharp  layers  would  be  expected  to  produce  at  long  ranges  are  not  likely  • 
to  be  observed  if  the  range  exceeds  the  layer  dimension.  Therefore,  in 
the  Chukchi  Sea  MIZ,  the  maximum  sound  level  fluctuation  may  level  off 
at  about  10  dB  when  the  range  is  increased  beyond  200  m. 

Another  method  of  predicting  the  acoustic  effect  of  thermal  layers 
was  derived  by  Potter  and  Murphy  (Ref.  4).  A  perturbation  calculation 
on  the  eikonal  equations  by  Bergmann  (Ref.  5)  was  extended  to  apply  to 
acoustic  transmission  in  a  medium  with  a  single  horizontal  layer  in  which 
the  deviation  of  the  index  of  refraction  from  unity  was  very  small.  Their 
results  applied  only  to  a  layer  with  a  temperature  variation  shaped  like 
the  Gaussian  error  curve.  The  effect  on  sound  level  predicted  by  this 
method  is  shown  in  Figure  68  for  three  locations  of  the  sound  source  with 
respect  to  the  layer.  Letting  o  be  the  standard  deviation  of  the  error 
curve  and  r  the  range,  the  predicted  sound  level  variation  is  given  by 

2 

L  =  -0.0109  AT  R  , 

cr 


where  R  is  the  value  given  on  the  curve  in  Figure  68,  and  AT  is  the 
magnitude  of  the  error  curve  in  degrees  Celsius. 


A  comparison  of  measured  sound  level  variations  with  those  predicted 
by  the  equation  above  shows  that  the  measured  variation  is  usually  much 
larger  than  the  predicted  one,  probably  because  temperature  layers  are 
not  usually  as  smoothly  shaped  as  the  Gaussian  error  curve.  The  depend¬ 
ence  on  the  square  of  the  range-to-a  ratio  should  be  noted  when  predicting 
fluctuations  for  ranges  longer  than  those  reported  here.  However  this 
increase  with  range  would  only  occur  if  the  layer  continued  throughout 
the  transmission  path. 


98  APL-UW  7311 


l 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


Figure  66.  Measured  and  predicted  sound  level  profiles  at  several  ranges 
source  at  SO  m. 
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Figure  68.  Sound  level  variations  predicted  by  Potter  and  Murphy. 


Statistical  Relationship 

When  many  small  thermal  layers  exist  it  is  not  possible  to  relate 
each  intensity  anomaly  to  a  certain  thermal  layer.  Instead,  we  sought  a 
statistical  relationship  between  the  average  size  of  the  thermal  layers 
and  the  resulting  average  intensity  variation. 

Using  the  data  taken  on  8  and  9  August,  when  temperature  profiles 
were  measured  at  4  stations  on  the  floe  every  15  minutes  for  6  hours,  the 
profiles  were  examined  to  determine  the  magnitude  and  horizontal  extent 
of  each  layer.  Layers  with  a  temperature  difference  of  less  than  0.03°C 
were  ignored  Figure  54  on  page  82  shows  the  distribution  of  the  layers, 
the  average  magnitude  of  their  temperature  difference,  and  their  average 
depth  dimension  for  each  day.  Because  of  the  large  spread  in  these 
values,  little  more  can  be  stated  than  that  the  layers  were  i  to  2  meters  ■ 
thick  with  a  temperature  difference  of  one-  to  two-tenths  of  a  degree 
Celsius. 

During  this  2-day  period  a  transmitter  was  placed  at  a  depth  of  33  m 
about  200  m  from  Camp  1  on  a  line  extending  beyond  Camp  2  sc  that  the 
transmission  passed  directly  under  Catip  2.  A  sound  intensity  profile  was 
measured  simultaneously  with  each  temperature  profile  at  Camp  1.  Each 
profile  consisted  of  sound- level  measurements  at  200  successive  depths  as 
the  receiving  transducer  was  lowered  or  raised.  The  pressure  amplitude 
wac  computed  for  each  reading,  and  the  average  and  the  standard  deviation 
of  the  amplitude  were  computed  for  each  profile. 
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atinnTnf  ”sult*«ejsh~n.in  Figure  69,  which  plots  the  standard  devi- 
tion  of  the  amplitude  divided  by  the  average  amplitude  for  the  profile 
The  figure  shows  a  greater  variation  on  9  August,  which  corresponds  to  ' 

the  layers  la/ei'S  ^  depth  dimension  of 


8 


Fzgure  69.  Statistical  variation  in  pressure  amplitude  for  each  profile. 

\ thermal  layers  varied  greatly  in  size  and  shape.  The  sound 
intensity  variation  depended  greatly  or.  the  proximity  of  the  acoustic 
source  to  a  large,  sharp  layer.  Therefore,  the  effect  of  the  layers  on 
acoustic  transmissions  can  only  be  stated  in  general  terms-  namely 
hat  layers  1  to  2  m  thick  with  a  temperature  difference  of  one-  tn  t-w 
for^nc^  3  d6-ree  Celsius  Produce  intensity  anomalies  averaging  +3  dB°" 
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VII.  VOLUME  REVERBERATION 


General 

A  knowledge  of  how  the  acoustic  scattering  strength  in  a  given 
region  of  the  sea  is  distributed  in  space  and  time  and  the  relationship 
between  the  scattering  strength  and  the  biological  scatterers  causing 
it  is  important  in  two  ways.  First,  the  more  direct  goal:  if  we  know 
the  acoustic  scattering  strength  in  a  given  region,  we  can  predict  to 
what  extent  active  acoustic  systems  such  as  military  sonars  or  fish¬ 
monitoring  equipment  may  be  limited  by  volume  reverberation.  Second, 
if  the  relationship  between  the  biological  parameters  and  the  measured 
scattering  strength  of  a  region  is  known,  it  should  be  possible  to 
monitor  the  biological  activity  in  the  region  with  relatively  simple 
acoustic  measurements,  and  to  build,  improve  and  test  biological  models 
of  the  region.  The  accomplishment  of  this  second,  more  long-range  goal 
would  enable  us  to  assess  more  accurately  the  effects  of  man's  activities 
on  the  local  ecosystem. 

The  direct  goal  is  clearly  attainable.  The  MIZ  studies  have  already 
contributed  sufficient  data  to  enable  reasonable  estimates  of  the  effects 
of  volume  reverberation  on  some  particular  acoustic  systems  in  a  given 
region  at  a  given  time,  viz.,  in  the  Marginal  Ice  Zone  in  summer  (Ref.  6). 
However,  the  present  work  is  only  a  beginning  toward  the  second  goal.  The 
relationship  between  the  acoustic  measurements  and  the  various  organisms 
obtained  in  the  biological  samples  is  now  being  studied.  We  are  also 
studying  biological  production  models  of  the  Arctic.  This  work  is  dis¬ 
cussed  further  in  Appendix  B,  a  paper  presented  to  the  1973  IEEE  Inter¬ 
national  Conference  on  Engineering  in  the  Ocean  Enivornment . 

The  first  measurements  in  the  ongoing  study  reported  here  were  made 
by  Moose  and  Shah  in  August  of  1971  (Ref.  1).  Measurements  of  the  acoustic 
volume  scattering  strength  were  made  at  35  locations  in  the  MIZ  at  105  kHz 
All  measurements  were  within  100  miles  of  Pt .  Barrow,  mostly  to  the  east 
in  the  Beaufort  Sea. 

In  1972  the  program  was  expanded  to  include  measurements  at  38  kHz 
as  well  as  105  kHz.  During  two  cruises  in  the  Chukchi  Sea  in  July  and 
August,  measurements  of  volume  scattering  strength  were  made  at  98  dif¬ 
ferent  locations  with  the  105  kHz  system.  Data  were  taken  with  the  38 
kHz  system  at  79  of  these  sites  and  at  two  other  locations.  CTD  measure¬ 
ments  were  made  at  most  of  the  sites,  and  at  35  locations  some  biological 
sampling  was  also  performed.  A  large  area  of  the  Chukchi  Sea  was  covered- 
as  far  north  as  73°00'N  at  169°14'W,  as  far  west  as  173°58'W  at  71°22»N 
(about  30  miles  east  of  Herald  Island),  and  as  far  south  as  66°50'N  at 
169  03 'W  (about  60  miles  north  of  the  Diomedes)  (see  Figure  22). 


104  APL-UW  7311 


I 

I 


T 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


seven ^  1972  measurements  were  taken  from  the  ship, 
seven  of  the  105-kHz  measurements  were  taken  from  helicopters--the  first 

:^Tmea?rmerS  S0  taken  in  thfe  Arctic  (or  “)"here  «l»e  to  our  knowl¬ 
edge).  These  data  were  taken  as  the  helicopter  hovered  about  8  to  10  ft 

iic°Hef°Pe!1>,WatK^’uleadS  °r  polynyas-  The  2ear  was  the  same  105-kHz  gear 

6  h  °28t  decS  1Pbrd  me*sulements  and  was  run  from  an  inverter  powfred 

ible  it  TiX,-n  6  heJlc°Pter-  The  technique  proved  quite  feas¬ 

ible.  The  use  of  helicopters  has  enabled  us  to  extend  our  range  of  on 

Sd  difficult regl°nS  °f  HeaVy  iCe  WhGre  Ship  °Perations  would  be  slow 


Experimental  Procedure 


-r 

J. 


I 

I 

I 

I 

I 
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71  anJhe-??PKrimental  method  has  been  described  previously  (Refs.  1  and 
7)  and  will  be  reviewed  only  briefly  here.  7  1  ° 

Inan  Both  the  105‘  and  38"kHz  systems  operated  in  essentially  the  same 

ab™tra  PU!Sef  t0.  °'6  ““  >0"S  —e  directed  downward  in 

bv  14”  fn!  ,K  t»Ct,TC  ”  IOr  the  105-kMz  systcm  and  a  beam  about  7" 

down  tn  38-kllz  system.  The  returning  echoes  were  amplified  beat 

pndV  ?  kH^’  and  recorded  on  magnetic  tape.  They  were  also  rectified 
and  displayed  on  paper  echogram  charts.  The  tape  recordings  were  then 
lgitized  and  analyzed  to  determine  the  average  acoustic  volume  scatter!™ 
strength  in  up  to  10  selected  depth  intervals  (Refs  8  and  91  Thp  d^nth8 
intervals  selected  for  each  measurement  site  were  "ho«e^har  the  ™hP 
charts  indicated  would  best  disclose  any  ll^  s^ctu^e  JHe  £££“ 

Thuf  “statistic']  flle  lnVOlVed  the  averaging  of  some  100-500  pings, 
tem^t-  1  err°rS  WerC  negUgible  (^f.  9)  compared  to  other  sys- 

Jertn  *tUnCerJairtleS'  tfe  estimate  the  uncertainty  in  the  volume  scat¬ 
tering  strength  of  any  given  layer  is  +3  dB  at  -50  to  -75  dB  (ref.  1  m"1). 

Figures  70  through  84  are  15  echogram  charts  accompanied  by  bar 
depth' interval  s'^Each'f^'  scatterln8  strength  in  the  corresponding 

£  variety'inlhelesults'obtainedl6''  “*  **“  ~P™t'"1f 

the  dltble  111  1iStS  thC  15  data  files  Presented,  the  location  at  which 
the  data  were  taken,  the  frequency,  the  lowest  depth  analyzed  the  maxi 
mum  volume  scattering  strength  in  any  layer  below  10  m  the  w^teS 
average  volume  scattering  strength  in  all  layers  below  10  m  and  the 
figure  displaying  the  data  file.  The  weighted  averlg^s,  is  given  by 


( 
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where  di  is  the  strength  in  decibels  (ref.  1  m'1)  in  the  ith  layer,  6^ 
is  the  depth  interval  of  the  ith  layer,  and  the  summations  are  taken 
over  all  layers  below  10  m. 

The  top  10  m  of  data  (the  top  11-12  m  of  the  ocean,  since  the  trans¬ 
ducer  was  generally  1-2  m  below  the  surface)  have  been  eliminated  from 
consideration  because  of  possible  interference  from  the  ship's  hull  or 
from  bubbles  caused  by  the  ship.  This  interference  was  most  noticeable 
in  the  38-kHz  data,  presumably  because  of  the  higher  sidelobes  in  the 
beam  pattern  of  the  38-kHz  transducer. 


Table  III 


Date 

(1972) 

Tape 

6 

File 

Lat.  N. 

Long  .  W . 

Freq. 

(kHz) 

Lowest 

Depth 

Analyzed 

(m) 

Maximum 

Sound  Scattering 
in  Any  Layer 
Below  10  m 
(dB  ref.  1  m"l) 

Average 

Sound  Scattering 
in  All  Layers 
Below  10  m 
(dB  ref.  1  m"l  ) 

Figure 

No. 

7/28 

5-2 

70*51 ' 

159°34 ' 

105 

28 

-71 

-77 

82 

j 

7/28 

5-3 

70°57' 

159°46 ' 

105 

53 

-77 

-80 

83 

( 

7/28 

5-4 

71°01 ' 

159°55 ' 

105 

52 

-80 
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81 

» 

7/28 

5-10 

71°3]  ' 

1 5  9  ^  C ' 
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51 

-75 
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80 

I 

7/28 

5-10a 

71°31 ' 

159°02 ' 

105 

51 
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80 

8/13 
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71°32' 

157*09* 

38 

110 
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-66 

76 

8/15 

6-19 

71°47' 

165°02 ' 

38 

40 

-58 

-60 

77 

1 

8/15 

6-20 

71°54 ' 

165*24 ' 

38 

40 

-60 

-63 

78 

8/13 

7-3 

71*29 ' 

157*02* 

105 

112 

-84 

-89 

72 

8/13 

7-4 

71°32 ' 
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105 
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-68 

73 

1 

8/13 

7-5 

71°37' 
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70 
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-73 

70 

| 
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7-6 

71°40' 
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63 
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71 

8/14 

7-12 
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i05 

62 

-66 
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84 

8/15 

7-16 

71°35' 
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ICo 

42 

-78 

-83 

79 

8/15 
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rr 

-61 
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74 
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mm 
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75 
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Figure  81.  Tape  5  File  4,  105  kHz ,  28  July  1972. 
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Figure  82.  Tape  5  File  2,  105  kHz ,  28  July  1972 
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Figure  83.  Tape  5  File  3 ,  105  kHz 3  28  July 
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Results 

Figures  70  and  71  are  typical  of  many  data  files  obtained  on  the 
cruises  and  help  to  give  aa  accurate  picture  of  the  average  overall  re¬ 
sults.  The  average  results,  however,  are  not  as  interesting  as  several 
specific  features  of  the  data  obtained  on  these  cruises. 

The  data  varied  greatly  over  very  short  distances.  A  good  example 
of  this  variation  can  be  seen  by  comparing  Figures  72  and  73.  These 
measurements  were  taken  only  3.8  n.  mi.  and  2  hr  apart,  with  little 
change  in  depth  or  ice  conditions  (about  6  oktas);  yet  there  was  a 
tremendous  increase,  about  21  dB,  in  acoustic  scattering  because  of  a 
dense,  shallow  layer  at  the  latter  site.  Clearly,  a  great  deal  of 
averaging  must  be  done  to  get  a  good  estimate  of  the  biological  activity 
in  an  area. 

data  taken  at  38  kHz  and  at  105  kHz  at  the  sane  location  showed 
practically  the  sane  scattering  strength.  (Compare  Figures  73,  74  and 
75  with  Figures  76,  77  and  78,  respectively,  which  were  taken  at  the 
same  sites.)  This  was  a  surprising  result.  Most  biological  scattering 
was  expected  to  come  from  organisms  such  as  plankton,  which  are  much 
smaller  than  the  acoustic  wavelengths  employed.  Therefore,  Rayleigh's 
law  should  have  applied,  and  the  acoustic  scattering  strength  at  38  kHz 
would  be  expected  to  be  approximately  18  dB  below  that  at  105  kHz.  There 
are  three  possible  explanations  for  this  phenomenon:  1)  the  scatterers 
were  large  compared  to  typical  plankton,  2)  they  were  closely  packed  so 
as  to  scatter  in  agglomerations  comparable  to  a  wavelength  (about  4  cm 
at  38  kHz),  or  3)  the  scattering  at  38  kHz  was  larger  than  normal  because 
bubbles  that  resonate  at  38  kHz  were  present  in  the  medium,  presumably 
attached  to  the  biological  species;  such  resonance  would  occur  at  a 
depth  of  25  m  for  a  bubble  about  0.05  cm  in  diameter. 

The  air  bubble  hypothesis  is  an  attractive  one  because  the  presence 
of  air  bubbles  varying  in  size  but  including  those  resonant  at  38  kHz 
would  tend  to  equalize  the  relative  effect  of  the  Rayleigh  scattering  at 
the  different  frequencies.  At  38  kHz,  one  38  kHz-resonant  bubble  per 
cubic  meter  would  produce  a  volume  scattering  strength  of  about  -31  dB 
(ref.  1  m-1).  The  possibility  of  scattering  from  gases  thus  would 
becloud  any  effort  to  determine  the  size  of  plankton  and  micronekton  bio¬ 
masses  by  comparing  their  acoustic  scattering  at  two  or  more  frequencies. 

In  general  the  38-  and  105-kHz  results  differed  less  than  the  com¬ 
bined  uncertainty  of  the  measurements.  The  scattering  strength  at  38  kHz 
was  not  always  as  large  as  that  at  105  kHz.  For  example,  the  definite 
scattering  layer  at  about  8  m  in  Figure  79,  which  was  taken  at  105  kHz, 
shows  a  volume  scattering  strength  of  -61  dB;  the  same  layer  had  a 
strength  of  -67  dB  at  38  kHz.  The  theoretical  Rayleigh  difference  of 
18  dB  was  never  observed. 
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It  should  be  noted  that  the  measurements  at  the  two  frequencies  were 
usually  made  within  a  few  minutes  of  each  other,  but  they  were  not  made 
simultaneously  because  of  mutual  interference  (reception  of  the  output 
pulse  of  the  other  system).  Figure  80  emphasizes  the  need  to  make  com¬ 
parative  measurements  on  a  specific  biomass  simultaneously.  The  dotted 
line  shows  the  average  acoustic  volume  scattering  strength  in  the  last 
206  pings  of  the  recording  when,  as  can  be  seen  on  the  echogram,  the  scat¬ 
tering  layer  at  about  34  m  concentrated  and  stabilized.  In  these  last  206 
pings,  the  average  scattering  strength  in  the  depth  interval  from  32  to 
36  m  jumped  4  dB  over  the  average  for  the  recording  as  a  whole.  There 
was  also  an  increase  of  5  dB  in  the  interval  from  20  to  32  m,  probably 
because  of  the  single  dark  object  at  about  23  m  near  the  right  edge  of 
the  echogram.  Meanwhile,  the  layer  near  15  m  became  less  pronounced  and 
yielded  about  5  dB  less.  In  the  future,  we  intend  to  synchronize  the  out¬ 
puts  of  the  two  systems  so  that  simultaneous  measurements  are  possible. 

Figure  81  is  of  particular  interest  because  of  a  sharp  echo  at  about 
8  m  on  the  echogram  that  corresponds  to  a  sharp  thermal  lager.  We  do  not 
know  whether  the  sharp  echo  was  caused  by  the  sudden  velocity  change  at 
the  thermal  interface  (that  is,  by  simple  partial  reflection  from  a  rho-c 
discontinuity)  or  by  organisms  that  congregated  at  the  interface.  A  rough 
calculation  indicates  that  the  physical  discontinuity  was  extreme  enough 
to  produce  the  observed  scattering.  Another  feature  of  Figure  81,  one 
typical  of  areas  where  there  were  warm  surface  intrusions,  is  the  congre¬ 
gation  of  biomass  in  the  warm  surface  water  and  the  relatively  little 
biomass  in  the  deeper  cold  water. 

Figure  82  (very  shallow  water ,  about  SO  m,  just  off  Pt.  Belcher) 
is  one  of  the  few  data  files  where  the  biomass  was  distributed  more  or 
less  evenly  throughout  the  volume.  This  relatively  homogeneous  dis¬ 
tribution  was  not  observed  farther  from  shore  and  was  probably  due  to 
in-shore  mixing  of  the  shallow  water.  Figure  82,  7  n.  mi.  farther  away 
from  Pt.  Belcher  in  water  about  55  m  deep,  shows  the  development  of 
definite  layering,  particularly  if  the  returns  that  appear  to  be  layers 
of  plankton  and  those  that  appear  to  be  single  targets  are  considered 
separately.  Another  4.6  n.  mi.  away  from  shore,  in  water  57  m  deep,  was 
the  extreme  surface  thermal  and  biolayer  previously  discussed  (Figure  81). 

Figure  84  is  unusual  in  that  it  shows  a  substantial  scattering  layer 
(-66  dB)  around  48  m.  Typically,  scattering  layers  above  -70  dB  were  net 
found  in  the  lower  half  of  the  volume.  The  large  temperature  discontinuity 
(the  top  12  m  are  7.5°C  warmer  than  the  rest  of  the  volume)  and  the  usual 
accompanying  concentration  of  biomass  near  the  surface  in  the  warm  water 
make  the  presence  of  this  deep  scattering  layer  all  the  more  unusual.  We 
do  not  know  its  composition. 
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Biological  Sampling 

the  ?,°  cruis”’  89  "«  hauls  net  trawls  were  made.  Of 

itV Pl  f'r!1  W6re  taken  With  a  3/4-m  diam-  ring  net  (216-u  mesh! 

(1000-y  mesh)  at  zlTsit^6"  3  ^  SqU3re  NI°  midwater  trawl  net 

^  ™d  “"bSeaftUe  ^  *° 

^5bW£&%33S£ 

was  the  only  net  available  that  could  be  closed  At  thP  lact  t! 

net  sites  the  net  was  jury-rigged  to  pivp  at  iMCt  th  1  St  t  °  ring 
the  dpnthc  cnmninj  l  I  8gea  zo  8ive  at  least  some  control  over 

depths  sampled,  but  this  arrangement  did  not  close  completely. 

The  ring  net  became  clogged  with  diatoms  whenever  used  All  in 

sst> 

to  errors  in  ?he  estimarlon  ef  f  5he>'  “eTe  ^significant  compared 

was  maintained  as  close  to  2  kn ^possible. DUrlng  the  toWS’  the  sPeed 

samples^ere  n^tatnTe^tnt  tlZln.  -‘—P-tion  plankton 

wire  angles  greater  than  65°  (over  80  m  of  paid-out  cable!  diffir.nJ 
was  also  experienced  in  tripping  the  net.  dl{flcultY 
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The  biological  sampling  was  performed  by  Mr.  Miles  J.  Furnas,  a  stu¬ 
dent  in  the  Department  of  Oceanography  at  the  University  of  Washington.* 

A  cursory  examination  of  the  samples  shows  a  large  number  of  various 
species,  and,  in  particular,  a  large  concentration  of  copepods  in  certain 
layers  (see  Table  I  in  Ref.  6).  A  few  individual  fish  were  found,  arctic 
cod  and  sculpin,  but  most  of  the  single  acoustic  targets  could  not  be 
caught  efficiently  by  the  method  used. 


Other  than  the  very  few  fish,  the  strongest  acoustic  scatterers 
found  in  the  1972  net  samples  were  euphausiids.  The  only  known  measure¬ 
ment  of  the  scattering  cross  section  of  an  individual  zooplankton  was 
done  on  euphausiids  by  Beamish  (Ref.  10).  Beamish  found  "hat  euphausiids 
typically  have  a  backscattering  cross  section  of  about  1.4  *  10“^  cm2  at 
102  kHz.  Thus  it  would  take  60  to  600  euphausiids  per  cubic  meter  to 
produce  an  acoustic  layer  with  a  volume  scattering  strength  of  -60  to 
-50  dB  (ref.  1  m"*)  at  105  kHz--if  the  scattering  were  caused  mainly  by 
these  organisms. 


Discussion  and  Future  Work 

A  comparison  of  the  105-kHz  results  with  those  of  1971  does  not  show 
any  differences  that  can  be  attributed  to  an  annual  change.  Some  of  thm 
very  dense,  relatively  deep  layers  observed  in  1971  were  missing  in  1972. 
However,  the  1971  data  were  taken,  on  the  average,  closer  to  shore  and 
under  less  ice  cover;  also,  a  considerable  part  of  the  1971  data  was 
taken  to  the  east  of  Barrow  in  the  Beaufort  Sea  instead  of  to  the  west  in 
the  Chukchi  Sea.  It  is  therefore  impossible  to  tell  whether  the  changes 
observed  in  the  scattering  were  caused  by  annual  differences  in  the 
climate  or  oceanographic  and  ice  conditions. 

Although  many  different  organisms  were  found  in  the  samples  collected 
in  1972,  information  as  to  the  volume  swept,  catching  efficiency,  and 
catching  depth  was  too  uncertain  to  warrant  any  attempt  to  use  the  data 
quantitatively.  More  accurate  sampling  methods  are  being  prepared  in  con¬ 
nection  with  the  biological  modeling  work  described  in  Appendix  B.  This 
improvement  should  complement  the  reverberation  work  by  enabling  us  to 
collect  samples  from  a  specific  acoustic  scattering  layer  more  efficiently. 
In  addition  to  improving  our  biological  sampling  capability,  we  plan  to 
add  a  250-kHz  system  to  collect  acoustic  data  on  the  smaller  zooplankton 
observed  in  the  samples. 

Because  of  the  great  local  variability  in  the  data,  any  quantitative 
conclusions  will  have  to  be  based  on  a  large  number  of  measurements. 


*Now  a  graduate  student  at  the  University  of  Rhode  Island. 
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There  were  179  data  files  taken  in  the  Chukchi  Sea  in  July  and  August 
1972.  Since  then,  August  1973  measurements  have  added  126  data  files. 
Analysis  of  the  data  taken  in  the  Chukchi  Sea  in  August  1973  should 
enable  us  to  make  an  appropriate  quantitative  comparison  of  the  scattering 
levals. 

More  data  will  be  obtained  in  July  and  September  1974  from  the  same 
region.  Statistical  analyses  of  these  data  will  be  made  and  correlated 
with  environmental  variables.  In  addition,  the  volume  reverberation  data 
taken  in  August  1971  in  the  Beaufort  Sea  and  in  March  1973  in  the  Bering 
Sea  will  provide  comparisons  with  other  regions.  As  more  data  on  the  re¬ 
gion  are  collected  and  analyzed,  we  will  be  able  to  gain  more  insight 
into  the  relationship  of  the  acoustic  scattering  strength  to  the  season, 
water  temperature,  and  ice  cover. 
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AUTOMATIC  PROFILER 

The  automatic  profiler  was  developed  as  a  prototype  unmanned  system 
for  gathering  conductivity-temperature-depth  data.  Basically  it  is  an 
automated  CTD  profiler  that  can  be  deployed  at  a  remote  site  and  left 
unattended  to  record  data  at  a  preprogrammed  rate. 

The  profiler  utilizes  the  same  types  of  sensors  as  the  CTD  system 
used  in  the  MIZ  and  records  data  in  the  same  format.  All  winch  operations, 
punch  control  and  drop  timing  functions  are  internally  controlled  by  the 
system.  The  duration  of  unmanned  profiler  operation  is  limited  primarily 
by  the  consumption  of  battery  power  and  the  capacity  of  the  recording 
system. 

System  Operation 

A  functional  description  of  the  automatic  profiler  is  shown  in  the 
block  diagram  in  Figure  Al.  The  operation  of  the  system  is  controlled  by 
the  master  control  unit,  which  can  program  profiling  at  a  rate  of  one  per 
hour  or  two  per  hour.  A  manual  override  control  also  allows  activation 
of  the  system  at  any  time.  The  control  unit  performs  all  data  processing 
and  tape  formating  operations  for  the  system  in  the  same  manner  as  in  the 
CTD  system  (see  Ref.  1). 


WINCH /SLIP  RING  ASSEMBLY 


MANUAL /AUTO 

DROP 

CODE 

CONTROL 

RATE 

SWITCHES 

Figure  Al.  Block  diagram  of  the  automatic  profiler. 
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When  the  control  unit  starts  a  drop  sequence,  power  is  automatically 
supplied  to  the  sensors  for  a  1-min  period  before  the  data  system  is  acti- 
vated.  During  this  time  a  2-ft  length  of  spacer  tape  is  punched  to  sepa¬ 
rate  the  new  data  record  from  the  previous  one.  After  the  initial  1-min 
delay,  data  samples  are  recorded  for  15  sec  with  the  probe  in  its  starting 
position.  This  data  is  used  to  provide  a  constant  starting  depth  refer¬ 
ence.  A  signal  from  the  master  control  unit  to  the  motor  control  unit 

then  applies  power  to  the  winch  motor  and  the  probe  is  lowered  at  a  rate 
of  1/3  m/sec. 

Lowering  is  continued  until  the  probe  touches  the  bottom  or  exceeds 
a  preset  depth  (80-m  maximum).  At  this  time,  data  sampling  is  stopped, 
the  winch  motor  is  reversed,  and  the  probe  is  raised  at  the  same  rate 
to  its  initial  position.  A  second  2-ft  length  of  spacer  tape  is  punched 

to  complete  the  record  and  all  functions  are  reset  to  await  the  next  drop 
command .  * 

System  Description 

The  winch  assembly  and  probe  are  shown  in  Figures  A2  and  A3,  respec- 
tiveiy  The  winch  is  mounted  on  one  leg  of  an  8-ft  tripod,  carries  100  m 
of  cable,  and  includes  a  module  with  a  48-volt  dc  motor,  gear  box,  drum 
and  slip  ring. 

switrhheaSH°neH^ntai?!-t?e  senso”f  the  Attorn  detector,  the  depth- limit 
switch,  and  a  data  multiplexer  and  its  control  electronics.  The  probe  has 

an  overall  length  of  50  in.,  a  maximum  outer  diameter  of  8.5  in.  and  an 
m-water  weight  of  26  lb. 

To  deploy  the  profiler,  a  hole  at  least  9  in.  in  diameter  must  be 
made  through  the  ice.  The  ice  thickness  must  be  sufficient  to  support 
the  totai  system  weight  of  600  lb  but  not  exceed  the  hole  drilling  capa- 
binty  The  tripod  and  winch  assembly  are  then  mounted  over  the  hole  and 
the  probe  is  dropped  through  it.  In  normal  operation  the  probe  is  left 
just  below  the  bottom  of  the  ice  to  reduce  the  possibility  of  damage  if 
the  probe  should  rise  automatically  when  the  hole  has  refrozen. 

The  batteries,  paper  punch  and  control  units  are  placed  on  a  sled  or 
paiiet  next  to  the  tripod  and  covered  with  a  canvas.  The  whole  system 
can  be  assembled  and  tested  in  less  than  1  hour  after  the  ice  hole  has 
been  cut.  The  deployed  profiler  is  shown  in  Figure  A4. 
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Figure  A2.  Automatic  profiler 
winch  assembly. 


Figure  AS.  Automatic  profiler 
CTD  probe. 
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he  primary  limitations  on  the  unmanned  operation  of  the  profiler  a 
the  capacity  of  the  recording  unit  and  the  battery  power  supplied  to  the 
system.  The  paper  tape  recorder  can  handle  1000  ft  of  tape.  The  take-up 
system  however,  is  only  90%  efficient.  At  the  profiler  drop  rate  of  1/3 
m/sec,  the  capacity  of  the  take-up  system  is  approximately  55  minutes  of 
profile  data,  or  20  drops  to  a  depth  of  50  m.  °* 


,_n.  Profller  1S  presently  supplied  by  two  48-V  packs  of 

to  frieS‘  ■  seParate  battery  pack  is  used  for  the  data  system 

isolate  it  from  noise  generated  by  the  dc  winch  motor.  When  the  pro- 

er,iS  n0t  ln  °Peration.  the  power  drain  from  the  punch  battery  is  less 
than  10  M.  The  »otor  battery  has  no  „o»er  loss  since  1?  Is  gS/out  oC 
the  system  until  the  drop  is  started.  During  the  drop,  the  average  bat- 
ery  drain  by  the  punch  is  2.5  A,  of  which  approximately  2  A  are  used  by 
the  recording  system.  The  motor  current  during  the  drop  varies  from  3  to 
4  A,  depending  on  the  ambient  temperature  and  the  drop  frequency.  A  150- 

0  4  A°h  fromnfhaPPr^XimKtely  0,1  A'h  °f  charSe  from  the  punch  battery  and 
'  tj®  m°t0r  battery-  At  an  ambient  temperature  of  approxi- 

:  L°a  •/nd  \rate  0f  me  50-  ,lr°r  the  "Otor  battery  "ill 

™rl  ?ha"lon?ops  reChargi"8  “  neCeSSar>'  lnd  the  "ett.^ 
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System  Improvements 

The  weakest  feature  of  the  profiler  system  proved  to  be  recording  the 
data  on  paper  tape.  Not  only  did  the  tape  take-up  and  supply  capacity 
limit  the  number  of  unattended  drops  but  the  punch  was  also  susceptible 
to  jamming  and  paper  tearing.  In  addition,  the  punch  consumed  relatively 
large  amounts  of  power. 

The  method  best  suited  for  data  recovery  of  this  type  would  appear 
to  be  data  telemetry.  This  method  would  eliminate  the  capacity  and  power 
supply  problems  of  an  on-site  recording  system  and,  in  addition,  would 
allow  transmission  of  data  relating  to  any  profiler  equipment  problems 
that  might  require  attention.  The  addition  of  a  radio  link  would  enable 
the  profiler  to  be  interrogated  from  the  recording  site  at  any  time. 

To  achieve  real-time  transmission,  the  telemetry  system  and  the  re¬ 
cording  system  would  both  have  to  be  capable  of  operating  at  a  data  rate 
of  250  bits/sec  or,  alternatively,  a  suitable  memory  system  could  be 
included  in  the  profiler  to  store  the  data  for  transmission  at  a  later 
time  and  the  data  telemetry  rates  could  be  adjusted  as  desired.  The 
latter  method  would  also  prove  useful  in  recovering  data  if  more  than 
one  profiler  were  being  controlled  by  the  data  recording  station. 
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EVALUATION  OF  A  TECHNIQUE  FOR  PREDICTING  VOLUME  REVERBERATION 
FROM  BIOLOGICAL  SCATTERERS  IN  THE  OCEAN 


ABSTRACT 

This  paper  describes  a  technique  for  predict¬ 
ing  the  scattering  strengths  of  biological  scat¬ 
tered  in  the  ocean.  An  attempt  was  made  to  pre¬ 
dict  the  seasonal  variations  in  the  biomass  and 
volume  reverberation  (at  various  frequencies)  in 
the  Arctic.  A  dynamic  biological  population  pro¬ 
duction  model  has  been  constructed  which  converts 
environmental  input  parameters  into  a  prediction 
of  the  density  and  site  distributions  of  the  bio¬ 
logical  population  as  a  function  of  space  and 
time.  This  output  is  then  converted  to  a  predic¬ 
tion  of  the  volume  scattering  strength  by  use  of 
theoretical  and  empirical  relationships  between 
organism  size  and  distribution  and  acoustic  scat¬ 
tering  phenomena.  We  have  used  a  simple  model  to 
describe  the  intertrophic  population  dynamics  of 
the  biological  scattered.  The  model  simulates 
relationships  between  the  nutrient  concentration 
and  the  population  density  of  the  various  organ¬ 
isms.  The  feasibility  of  monitoring  the  biomass 
by  acoustic  techniques  coupled  with  computer 
simulation  is  also  discussed. 


INTRODUCTION 

The  major  source  of  volume  reverberation  in 
the  oceans  has  been  established  as  biological  or¬ 
ganisms  [1],  Non-biological  sources  such  as  dust 
and  sand  particles,  thermal  gradients,  and  turbu¬ 
lence  of  natural  or  man-made  origin  are  usually 
insignificant  contributors  to  the  volume  scatter¬ 
ing  strength  observed  at  sea. 

Many  studies  have  been  done  to  identify  the 
biological  organisms  responsible  for  sound  scat¬ 
tering,  as  well  as  to  determine  the  scattering 
properties  of  an  individual  organism.  Love  [2] 
and  Cushing  and  Richardson  [3]  have  made  measure¬ 
ments  of  the  target  strength  of  fish.  Love  has 
determined  a  description  of  the  target  strength 
of  an  individual  fish  as  a  function  of  its  length 
and  the  acoustic  wavelength. 

During  the  last  decade,  much  interest  has 
been  shown  in  the  acoustic  assessments  of  fish 
stocks  [4,5] .  The  recent  developments  in  high¬ 
speed  processing  of  sonar  information  render  the 
task  of  measurement  and  estimation  of  the  volume 


reverberation  easier.  However,  we  are  far  from 
knowing  the  volume  scattering  strength  in  differ¬ 
ent  oceans,  at  various  seasons,  depths  and 
frequencies. 

Figure  1  is  a  simplified  block  diagram  of  the 
dynamics  of  the  marine  ecosystem.  The  principal 
biological  components  are  phytoplankton,  and  her¬ 
bivorous  and  carnivorous  zooplankton,  with  nutri¬ 
ents  and  visible  radiation  as  the  driving  forces. 


Figure  1.  Marine  Ecosystem. 

There  is  interaction  among  the  biological  compo¬ 
nents,  which  collectively  form  biomass  to  produce 
volume  reverberation.  Successful  prediction  of 
the  bioacoustic  environment  can  lead  to  confident 
prediction  of  mean  volume  scattering  strengths. 

The  volume  scattering  strengths  are  estimated  from 
the  predicted  distribution  of  biological  scatter¬ 
ed  ;  therefore,  it  is  necessary  to  know  the  aver¬ 
age  body  size,  mass  and  the  number  of  different 
scatterers.  Based  on  the  body  size,  the  target 
strength  for  various  scatterers  can  be  estimated 
at  different  frequencies. 

Since  1971,  the  Applied  Physics  Laboratory  of 
the  University  of  Washington  has  been  actively  in¬ 
volved  in  the  measurement  of  volume  reverberation 
in  the  Arctic  as  part  of  the  Marginal  Ice  Zone 
Pacific  Project  (sponsored  by  the  Arctic  Submarine 
Laboratory,  Naval  Undersea  Center,  San  Diego)  [6]. 
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In  addition  to  Baking  quantitative  oeasurenents 
of  volune  scattering  strength,  biological  samples 
were  taken  to  identify  the  type  and  relative  num¬ 
bers  of  the  scatterers.  As  a  result  of  the  mea¬ 
surements  taken  in  the  Arctic  (Chukchi  Sea  and 
Beaufort  Sea),  and  because  of  the  apparent  sim¬ 
plicity  in  biological  annual  cycles,  the  lower 
Arctic  in  the  vicinity  of  Barrow  was  chosen  as 
the  initial  area  of  investigation. 

PRELIMINARY  POPULATION  MODEL 

Plankton  population  dynamics  is  influenced 
by  many  complex  biological  and  physical  phenomena. 
It  is  necessary,  therefore,  to  idealize  and  simp¬ 
lify  the  conceptual  model  so  that  the  result  can 
be  manageably  simulated.  Referring  to  Figure  1, 
the  state  variables  of  the  system  are  concentra¬ 
tions  of: 

phytoplankton  (p) 
herbivorous  zooplankton  (h) 
carnivorous  zooplankton  (c) 
nutrient  (n) 

The  visible  radiation  intensity  (with  temporal 
variations)  is  the  forcing  function. 

The  model  simulates  relationships  between  the 
state  variables  in  a  fixed  volume  of  water  (for  a 
homogeneous  area).  The  simulation  consists  of 
the  solution  of  a  system  of  coupled  partial  dif¬ 
ferential  equations,  one  for  each  of  the  state 
variables. 

Our  simulation  of  the  model  is  based  on  the 
following  assumptions: 

1)  The  photosynthesis  rate  is  proportional 
to  the  light  intensity  at  low  intensity, 
but  is  limited  when  the  saturation  rate 
is  reached  according  to  Steele's  [7] 
expression 

Ph.Ph(max)(I/Iopt)exp(l-I/Iopt), 

where  Pfj  is  the  photosynthesis  rate,  I  is 
the  incident  light  intensity  and  I0pt  is 
tLe  light  intensity  which  gives  maximum 
photosynthesis  rate  Pj1(max) . 

The  light  intensity  I  at  depth  z  as  a 
function  of  incident  surface  radiation  I0 
may  be  expressed  by  the  Beer-Lambert  law 

as 

I  -  IQexp(-kz) 

where  k,  the  coefficient  of  extinction, 
is  assumed  constant.  Self-shading  effects 
may  be  included  by  making  k  a  function  of 
phytoplankton  concentration;  that  is, 
extreme  increases  in  phytoplankton  would 
inhibit  tl  "  sunlight  and  reduce  photosyn¬ 
thesis  production. 

2)  Inorganic  nitrogen  is  considered  to  be 
the  limiting  nutrient  affecting  the 


B2  APL-UW  7311 


photosynthetic  assimilation  of  carbon  [8], 
The  ratio  of  photosynthetic  assimilation  of 
carbon  to  nitrogen  is  assumed  to  be  6  [9], 

3)  The  Ivlev  [10]  predator-prey  relationship 
describes  the  grazing  and  predation  rates. 
The  Ivlev  expression  for  grazing  is 

W  ‘  "m.xI|-e‘kP]> 

where  M,nax  is  the  maximum  grazing  rate,  p 
is  the  phytoplankton  (prey)  concentration 
and  k  is  a  constant. 

4)  The  fractional  loss  of  phytoplankton 
biomass  due  to  sinking  is  constant. 

5)  Respiration,  excretion  and  mortality 
rates  are  constant  for  each  trophic  level. 

6)  A  euphotic  zone  exists  from  the  surface 
down  to  20  meters  (based  on  50%  light 
level  at  10  meters),  with  homogeneous 
distribution  of  the  bio-concentrations 
in  this  zone. 

7)  Water  temperature  is  constant  at  0*C. 

8)  Assimilation  efficiency  of  herbivores 
is  60%  [11],  and  for  carnivores  is  85%. 

9)  Nutrient  flux  is  described  by  the  Fickian 
expression  for  diffusion, 

flux  »  -kz  •  dn/dz, 

where  kz  is  a  vertical  eddy  diffusivity 
constant,  and  dn/dz  is  the  nutrient 
vertical  gradient  [12], 

The  coupled  partial  differential  equations 
used  in  the  simulation  are: 

3p/3t  -  p(Ph  -  RP  -  SINK)  -  Wj  •  h 

3h/3t  »  h(0.6Wj  -  RH  -  MOLT)  -  *2  •  C 

3c/3t  »  c(0.8SW2-  RC  -  MORT) 

3n/3t  -  h  •  EXC1  ♦  c  •  EXC2  -  kz  •  dn/dz  -  (p  •  Ph)q 

The  symbols  are: 

RP  -  respiration  rate  of  phytoplankton 

RH  -  respiration  rate  of  herbivores 

RC  -  respiration  rate  of  carnivores 

SINK  -  sinking  rate  of  phytoplankton 

MOLT  -  mortality  rate  of  carnivores 

*1  *  JTazing  rate  of  phytoplankton  by  herbivores 

1*2  -  predation  rate  of  herbivores  by  carnivores 

EXC1  -  excretion  rate  of  ammonia  by  herbivores 

EXC2  -  excretion  rate  of  carnivores 
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q  -  conversion  factor  between  the  Mount  of  carbon 
assimilated  and  equivalent  amount  of  inorganic 
nitrogen  consumed. 

The  simulation  of  this  simple  model  consists 
of  solving  the  four  partial  differential  equations, 
by  the  elementary  forward  time-step  method.  That 
is,  if  we  represent  the  four  equations  by 

3F.it) 

— |t“"  fi  i  - 

(where  the  F^  are  the  four  concentrations  and  $(t) 
represents  the  driving  functions);  then, 

3F,(t) 

Fj(t*At)  ■  Af- ♦  F^t);  i  -  1 ....  ,4. 

In  this  way,  given  the  initial  concentrations  of 
inorganic  nitrogen,  phytoplankton,  herbivorous 
and  carnivorous  zooplankton  and  their  functional 
interdependence,  we  can  project  their  concentra¬ 
tions  forward  in  time. 

The  stability  and  accuracy  of  the  method  are 
preserved  if  the  changes  in  the  variables  during 
time  At  are  small  enough  so  that  the  first  order 
term  in  Taylor's  expansion  is  adequate.  This  ele¬ 
mentary  method  of  solution  has  proven  sufficient. 

PRELIMINARY  ACOUSTIC  MODEL 

The  acoustic  volume  scattering  strengths  are 
obtained  using  the  biological  scatterer  density 
predicted  by  the  population  model  simulation.  It 
is  therefore  very  important  to  know  the  target 
strengths  of  individual  scatterers  as  a  function 
of  frequency.  The  target  strength  of  a  distribu¬ 
tion  of  organisms  which  scatter  incoherently 
relative  to  one  another  is  given  by  [1] 

T  -  10  log  C  N(S)o(S,X)dS, 

0 

where  o(S,X)  is  the  average  scattering  cross  sec¬ 
tion  at  acoustic  wavelength  X  and  N(S)dS  is  the 
number  of  organisms  per  unit  volume  whose  length 
is  between  S  and  S*dS. 

Computation  of  the  Target  Strengths 

Cooney's  [13]  expression  was  used  for  the 
small  targets  (<1  cm  in  length)  and  Love's  ex¬ 
pression  was  used  for  fish. 

Cooney's  expression  for  the  target  strength 
of  very  small  targets  based  on  Rayleigh  scatter¬ 
ing  for  the  individual  plankton  is 

if2j2  3  2 

TS  •  10  log  (1  ♦  yp)  , 

where  TS  is  the  target  strength  in  decibels,  T  is 
the  volume  of  the  "ideal"  animal,  X  is  the  wave¬ 
length  of  the  incident  sound,  and  w  is  the  cosine 
of  the  angle  between  the  scattering  direction  and 
the  reverse  direction  of  the  incident  wave;  |i>l 
for  backscattering. 


Love  [2]  has  developed  an  empirical  equation 
that  approximates  the  dorsal-aspect  target 
strength  (Tp)  of  an  individual  fish  in  the  L/X 
range  of  interest  for  most  sonar  application, 
where  L  is  the  length  of  fish,  and  X  is  the 
acoustic  wavelength. 

Love's  expression  for  the  dorsal-aspect 
target  strength  of  an  individual  fish  is 

TD  •  19.1  log  L  ♦  0.9  log  X  -  23.88, 

where  L  and  X  are  in  meters,  and  Tp  is  in  deci¬ 
bels.  The  equation  is  a  good  approximation  to  the 
dorsal -aspect  target  strength  of  an  individual 
fish  for  0.7  <,  L/X  <  90. 

SIMULATION  OF  THE  MODEL  AND  THE  RESULTS 

The  model  was  simulated  to  predict  the 
concentrations  of  phytoplankton,  herbivorous  and 
carnivorous  zooplankton  (in  milligrams  carbon/ 
meter^)  and  inorganic  nitrogen  (in  pgAt  N/liter) 
from  spring  to  fall  on  a  daily  basis.  The 
volume  reverberation  is  predicted  at  a  frequency 
of  100  kHz,  the  frequency  used  for  the  measure¬ 
ment  of  volume  reverberation  in  the  Arctic  in 
1971  and  1972. 

Results 

Figure  2  shows  some  of  the  results  for  a 
typical  set  of  inputs  as  shown  in  Table  I.  There 
are  early  spring  phytoplankton  blooms  which  are 
limited  by  the  inorganic  nitrogen  concentrations. 
The  blooms  are  depleted  by  the  herbivorous  zoo¬ 
plankton  stock,  which  in  turn  is  consumed  by  the 
carnivores.  The  typical  predator -prey  relation¬ 
ship  is  exemplified  by  the  oscillatory  behavior 
of  the  plankton  variables.  The  volume  scatter¬ 
ing  strength  varies  between  -85  dB  and  -70  dB. 

The  peaks  in  scattering  strength  coincide  with 
peaks  in  carnivore  concentration.  The  scatter¬ 
ing  data  at  105  kHz  obtained  experimentally  in 
1971  are  also  shown.  An  elementary  parameter 
sensitivity  analysis  resulted  in  the  following 
observations : 

1)  A  10%  variation  in  inorganic  nitrogen 
concentration  does  not  have  a  signif¬ 
icant  effect  on  the  scattering  strength. 

2)  Increasing  the  respiration  rates  of  the 
plankton  by  10%  also  affects  the  scater- 
ing  insignificantly. 

3)  Increasing  the  grazing  or  feeding  rates 
by  10%  increases  the  scattering  by 
about  2  dB  at  the  peak. 

4)  Doubling  the  initial  concentration  of 
the  state  variables  resulted  in  an 
increase  in  scattering  of  about  5  dB 
during  the  first  4  weeks,  after  which 
the  model  stabilized  with  only  slight 
variations  in  scattering. 
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TabU  I. 

Initial  Concentration  Quantity  Unit 

phytoplankton  10  mg  carbon/a3 

herbivores  5  ag  carbon/a3 

carnivores  3  ag  carbon/a^ 

inorganic  nitrogen  3  yg  At  N/liter 

Respiration  Rates  (fraction  per  day)  [14] : 

phytoplankton  0.15 

herbivores  0.13 

:amivores  0.05 


Figure  2. 


DISCUSSION 

The  siaulation  of  this  siaple  aodel  has  shown 
that  with  the  knowledge  of  certain  biological 
paraaeters,  voluae  reverberation  can  be  predicted 
with  confidence.  Accuracy  of  the  aodel  is  depen¬ 
dent  upon  a)  the  validity  of  the  aodel  and  b)  the 
accuracy  in  the  knowledge  of  the  biological  rate 
paraaeters . 

Though  the  acoustic  scattering  strengths  ob¬ 
tained  froa  the  siaulation  fall  within  the  range 
of  the  experiaental  data,  this  aust  be  viewed 
critically,  and  not  as  an  indication  of  the  suc¬ 
cess  of  the  aodel.  Experiaents  should  be  designed 
to  obtain  the  critical  data  needed  for  operation 
and  verificat!  n  of  this  aodel  and  its  simulation. 
This  would  entail  biological  saapling  and  acoustic 
and  physical  aeasureaent  procedures  to  be  per- 
foraed  over  sufficient  tiae  and  averaged  over 
sufficient  area  to  be  a  critical  test  of  the  bio¬ 
logical  mode,  its  siaulation,  and  its  acoustic 
consequences . 

The  current  interest  in  modeling  the  oceanic 
upwelling  regions  of  the  world  has  resulted  in 
large-scale  simulation  programs  under  the  Inter¬ 
national  Biological  Program.  These  involve  a 
team  of  specialists  simulating  and  verifying 
models  on-site.  These  long-term  programs  should 
provide  better  understanding  of  the  population 
dynamics  and  more  accurate  rate  parameters. 

Studies  should  be  carried  out  which  entail  a  bio¬ 
logical  and  acoustic  investigation  in  a  small 
oceanic  area,  conducted  over  a  period  of  time,  to 
obtain  a  better  understanding  of  the  relationships 
between  biological,  acoustic  and  physical  oceano¬ 
graphic  variables.  The  knowledge  and  experience 
gathered  from  investigations  in  one  area  should 
lead  to  prediction  models  for  other  areas  of 
interest. 


Maximum  photosynthesis  rate  >  1.00  per  day 

I  t  *  0.5  gm-cal/cm2/min 

Extinction  coefficient  of  sea  water  >  0.05/m 

Sinking  rate  >0.05  per  day 

Mortality  rate  of  carnivores  ■  0.10  per  day 

Molting  rate  of  herbivores  >  0.05  per  day 

Maximum  grazing  rate  >1.5  per  day(») 

Maximum  predation  rate  >0.6  per  dayM 
Ivlev  constant  for  herbivores  -  0.015^*^ 
Ivlev  constant  for  carnivores  >  0.029^“) 

Individual  [14] 


Average 

Average 

Target 

length 

mass 

strength 

(-) 

(mg  carbon) 

(dB) 

Herbivores 

5 

0.039 

-110 

Carnivores 

15 

1.0 

-85 

Light  data  were  obtained  from  English's 
measurements  in  the  Arctic  [15]. 

yy 

*'  Unpublished  data,  M.  Jawed,  Dept,  of 
Oceanography,  University  of  Washington 
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